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NON DAECY PLOW - A PEESPECTIVE 
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IHTROrOCTION 


Ron Darcy Plow and its Importance 

The flow of water which follows the famous 
Darcy's law (1856) expressed as 

V = lei . . . . . ( 1 ) 

v/here v = Superficial flow velocity 
i = Plydraulic Gradient 
k = Constant of Proportionality called 
coefficient of permeability. 

is called the Darcian flow, while the flow not expressible 
by above eqn, (1) can be designated as non Darcy flow, 

Darcy's law has made a tremendous impact on the 
understanding of the movement of water in solids and 
porous media. Darcy's law when combined with the equation 
of continuity leads to heat conduction or diffusion type 
equations which are mathematically solvable for various 
boundaiy conditions. Any refinement in the Darcy's 
equation either may result in an unsolv.able nonlinear 
differential equation or make the boundaiy conditions 
difficult for mathematical treatment . But the mathe- 
matical difficulties have not deterred a study of non 
Darcy flow behaviour. 



It is to be noted here that Darcy (l856) established his 
experimental law under the following six conditions, 

1. Mediijm size sand 

2» The sand bed was stable; that is there was no 

structural or particle orientation changes during 
flow, and movement of fines was negligible, 

3. The sands tested were of inert type i.e», there was 
no surface force to affect the pemeant, 

4. The flow was under isothemal condition 

5. The testing was conducted over medim range of 
gradients 

6. The flow Was under steady state condition. 

If one’s testing conditions are different from Darcy’s, 
deviations from Darcy's law might be expected. 

Consequences of non Darcy behaviour are of potential 
interest in several disciplines. Drainage and water movement 
in clay soils are of primary importance to soil sciences, 
hydrology and Soil Engineering, Consolidation of clay depoeiiti? 
caused by water migration is of utmost interest and ftindamenibel 
importance to Soil Mechanics, The consideration of non 
Darcy flow during consolidation may perhaps remove many 
descripencies. Between observed and predicted results based 
on present Darcian theory. The existence of threshold 
gradient has also important implications in consolidation of 
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saturated clays and for water sealing properties of 
clays (Swartz endru be r 1968). Non Darcy behaviour 
will also drastically affect the flow through unsaturated 
soils which is of primary importance to agriculture and 
hydrology. 

Some of the applications of non Darcy flow in 
practical problems which have been solved by different 
authors may be mentioned here, Hansbo (I960), 

Elnagger et al (1971) and Parlange (1973) dealt with 
the consolidation of soils depicting non Darcy flow 
behaviour. Bondarenko (1968) calculated the drainage 
spacing for soils showing the phenomenon of initial 
gradient. Swartz endruber (1968) gave the head, 
distribution in an acquifier having non Darcian behaviour. 
Slepioka (1961a) studied the discharge from wells 
placed in a non D&rcian aquifer. 

Non Darcy flow of v/ater through soils is generally 
observed under two circumstances and they are 

(1) Blow at high gradients - the main cause 
of deviations being inertial effects. 

Blow at low gradients - the main cause of 
deviations here is due to the predominance 
of soil-water interaction at the molecular 
level and the associated effects. These 
are explained later in detail. 


( 2 ) 
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Darcy's law is generally observed to be perfectly 
valid for intermediate range of gradients . In this 
report, the non Darcy flow which is observed only at 
the low gradients, is discussed unless mentioned 
otherwise. 
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Chapter 1 

'■ CAUSES OE UQU-DARCIAJT BEHAVIOUR 

There are several hj^'potheses used to explain 
the non Darcian flow in the laminar region and these 
can be grouped under the following six heads 

(1) Hon Uewtonian liquid viscosity hypothesis 

(2) Bingham Plastic fluid and '’Quasi-Crysta- 
-lline" water structure hypothesis 

(3) Electro kinetic coupling hypothesis 

(4) Wall force hypothesis 

(5) Domain structure hypothesis 

(6) Transient particle arrangement hypothesis 
A brief description of each of the above 

hypotheses will be given below, 

(1) Hon Newtonian liquid viscosity hypothesis ; 

Under this hypothesis non Darcian flow is 
supposed to originate from non Newtonian liquid 
viscosity caused by modified water structure induced 
by clay water interaction, low (1961) first developed 
this hypothesis and later it was extensively used by 
different authors to explain their experimental data. 
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(Swartzendruber, 1962a, 1962b, 1966, 1967? Kutilek 

1964, 1965, 1967a, 1967b, 1969 and 1971? Kutilek 

and Salingerova 1966, 1967? Karadi & Nagy 1961). 

Von Engelhardt and Tunn (1955), Lutz and Kemper (1959), 

Hansbo (I960), Hadas (1964^ and Thames (1966) also 

explained their experimentally observed non Darcian 

behaviour on the basis of increased viscosity due to, 

surface activity of soil minerals. Kutiliek (1964 & 

1967 ) proposed that some part of water in a clay -water 

system behaves as a Newtonian fluid whereas the other 

paft behaves as non-Newtonian liquid. A detailed 

discussion about the changed viscosity and density 

of water in a clay-water system is given later in 

of 

this report. Experimental results of each^the above 
mentioned researchers are also given in a separate 
chapter, 

(2) Bingham Plastic fluid and "Quasi crystalline" 
water structure hypothesis 

Bondarenko (1968), Churyov and Gorokhov (1970), 
Nerpin and Beryagin (1969), Millar and Low (1963) and 
many others believe that water in capillaries has the 
properties of Bingham body. ITerp::,n or Deny, agin and 
Bondarenko measured the yield shear stress of capillary 
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water of the order of 10 to 10 dynes/cm , and 
hence, according to them initiation of flow in narrow 
capillaries of soils needs certain hydraulic gradient 
and thereby violating Darcy's Law at low gradients. 
Bondarenko and others attributed this property of 
capillary water to the hydrogen bonding and quasi- 
crystalline water structure in the capillaries. 

(3) Electrokinetic Coupling Hypothesis 

Electrokinetic coupling is the interaction that 
lAccurs between viscous and electrical flows in 
membrane materials such as clays (Abramson, 1936, 
Overbreek 1952). When a liquid is forced thro\;igh 
clay, by a hydraulic gradient, the coupling gives rise 
to an induced streaming potential that causes an 
electro-osmotic flow opposite to the flow caused by 
hydraulic gradient (Bwartzendruber 1966, 196?) resulting 
thereby in a reduction in the flow i^ate as predicted 
from any theory based on Darcy's law (e.g. Poiseulle's 
law or Kozeny Carman's equeiion) Swartzendruber (1967) 
gave the velocity -gradient relationship as 

V = ki - Gig . . . . . (2) 

or V = (k - 0 EA)i . . . . (3) 
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where ± =. :=■ = hydraulic gradient 
Jj 

i^ = ^ = streaning potential gradient, 

E is the streaming potential 
across the sample of length L. 

If E/h is assumed as constant (i..'.-;:r'mison 1934, 
Overbreek 1952 ), v is still found to be directly 
proportional to hydraulic gradient i,e Darcian flow 
is observed but Bull and Gortner (1932) and Swartzen- 
druber (1966) found that E/h decreases with ’h’., and 
hence, with hydraulic gradient, thereby making the flow 
non Darcy. Henniker (1952) also observed a retar- 
dation in flow at low gradients and he also used 
streaming potential to explain it. Olsen (I960) 
calculated the effect of electrokinetic coupling on 
the flow rate in percent and he found that the coupling 
influence is negligibly small (less than 1^), 

(4) Wall force hypothesis; 

Taking hints from Martin (I960), IClousner and 
Kraft (1965a, 1965b, 1966a and 1966b) proposed a 
capillary model with a force E normal to the surface 
(resulting from the electi’ical field from the crystalline 
structure of soil minerals) and hypothesised that this 
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force F brings a shear force which increases as the 
solid particles', is approached. With some theoretical 
deduction they were able to show that the existence of 
the wall force can give rise to the non linear velocity 
gradient relationship. For their analysis Klausner 
and Kraft assumed that there is no variation of 
viscosity and density of pore fluid. So the entire 
non Darcy effect is assimied to be due to existence of 
wall forces alone. 

(5) Domain Structure Hypothesis; 

Michaels (1959) and Quirk (1959) suggested;.'" 
that the primary particles in a clay mass may be 
arranged in groups such as aggregates, packets or 
domains, and that the total porosity may be distributed 
among inter and intra group components. If such 
grouping arrangement of particles exists in clays, 
the flow channels surrounding partic..uvv groups probably 
will be considerably larger than those- passing through 
the groups and betvreen the individual particles. If 
this sort of domain or cluster aotructure is accepted 
then the Kozeny-Carman' s permeability equation based 
on parallel capillary tube model and Darcy's law will 
be cdmplet’ely invalid. The failure o:^ Kozeny -Carman' s 
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equation might mean the failure of parallel capillary 
tube model and might warrant a better model. 

Following the domain structure or cluster 
concept Olsen (1961) derived a flow rate equation 
which reads as 

!S2S ^^ 2/2 h-UEl (4) 

%ozeny 4/3 

(1 + 


where = Flow rate for soil with domain 


or cluster structure 


^ozeny 




KqS m 1 + e 


i A 


• • • • (5) 


where e = void ratio, = shape factor, 

S = specific surface, Y and 
are viscosity and density of liquid, 
and A is area of flow. 

Parameters in equation (4) ar^. explained in 
Figs, 1 and 2. 


Using eqn. (4) Olsen calculat.::d the values of 


discrepency ratio (• 


"NFS 


-) for different porosities 


'%Cozen,Y 

with a definite cluster geametry (Fig. 3). 


Compared to Kozeny-Carman’ s equation Olsen's 












12 


equation (eqn. 4) may give better result for flocc- 
ulated and aggregated soils and might explain some of 
the "apparent" deviations from Darcy's law, 

(6) Transient Particle Arrangement Hypothesis : 

Michaels and lin (1954), Martin (I960), Olsen 
(I960, 1961), Swartzendruber (1962b), Miller and Low 
( 1963 ), Zaslavski ( 1964 ), Blackmore and Marshall 
( 1965 ), Mitchell and Younger (1967), Marshall (1968), 
Gairon (1968), Miller, Overman and Perverly (1969), 
Kutilek ( 1969 ) and Movak (1970) suggest that the 
particle re -orient at ion or the change of the geometric 
arrangement of particles in the soil sample (i.e., 
fabric or structural change) is the most important 
effect in non Darcy flow, 

Michaels and lin's (1954-) data (quoted in Martin 
I960) demonstrate that clay fabric alone can alter 
the permeability at a given void ratio by 1000 percent. 
It is possible that the application of hydraulic 
gradient to a structure or fabric sensitive clay will 
change its fabric appreciably (void ratio may or may 
not remain same during this process) giving rise to 
an observed non Darcy behaviour Martin (I960) explained 
Von Engelhard t and Tunn's (1955) non Darcy data in terns 
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of reversible void plugging or movements of fines from 
one pairt of the sample to the other parts and 
concluded that mejor cause of non larsian behaviour 
is the aggregation or dispersion of fines. 

/ Olsen (I960, 1961) derived a relation between 
flow rates and orientation of particles in a soil mass 
and his equation -reads as 


%DPO 
^Ko 


zeny 


1 + 


1 


T1 4 - 4.-i£ - .1 
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■ ■( 6 ) 


Sin^© 


where <lj;]-ppQ = Plow rate for a particular orientation 
of particles 

'^Kozeny “ Plow rate' given by Kozeny-Carman’ s equation, 
Other parameters of the equation (6) are explained 
in Pig. 4. '■ _ 

The discrepancy., ratio ■ ■) is calculated for 


different particle orient‘at^cnt-and is shown in Pig. 


5. As seen-.:f-rom this at- the same porosity or 

void ratio, as the particle orientation changes (due 
to the application of hydraulic gradient) the flow 
velocity chaiges giving rise to an apparent non Darcy . 



Direct tOQ of Hydraalic 
Gradieit 


14 



I 





behavioux. 
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Swart zendruber (1962a and 1962b) hinted that if 
with increasing gradient the particles orient 
themselves in such a wa,y that effective cross section 
of flow area changes then one ¥70uld orpect a non Darcy 
behaviour. 

Miller and low (1963) considered the 
rearrangement of particles during their flow experiments, 
to be parallel to the flow path and the velocity- 
gradient relation was no more linear as the tortuosity 
deorsaafi-s* 

Zafilavski (1964) foimd that shift of particles 
cause hydraulic conductivity (slope of v-i curve) to 
be time and gradient dependent, 

Blackmore & Marshall (1965) .have shown that in 
a bentonite suspension the shift of particles causes 
a gradual decrease in spacing, the change being 
continuous and dependent upon applied hydraulic 
gradient and velocity. 

Mitchell and Yo\mger (4-967) observed that the 
behaviour of compacted spepimens was such as to 
indicate that particle migration as the more likely 
cause for non Darcy flow. 
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Marshall (1968) found that, when water flows 
through a bed of unstable material, the hydraulic 
conductivity and hence velocity gradient relationship 
can be affected by detachment of particles or by 
change of the volume of the material. 

Gairon ( 1968 ) , (quoted in Kutilek 1969 ) conducted 
experiments with bentonite clay pastes and according 
to his results the flow was less than proportional, 
at least at low concentration and be explains these 
deviations from the proportional Darcian flow mainly 
by the displacement of particles within the column 
of clay paste. 

Miller, Overman and Peverly ( 1969 ) observed 
that not only simple movement but also bending and 
flexing of particles or breaking of edge to surface 
bonds, to permit particle reorientation along the flow 
path may cause non Darcy flow. 

Kutilek (1969) analytically demonstrated the 
dependence of hydraulic conductivity and hence 
velocity-gradient relationship, upon the hydraulic 
gradient in suspension of platy particles. 

Novak (1970) explained his nonlinear- velocity 
gradient' relation for bentonite in terms of shifting of 
the position of bentonite particles conforming to the 
velocity and direction of flow. 



CHAPTER 2 


n ' 


BOUHID WATER AKI) ITS ROIE IH HOF FARCY BEHAVIOUR; 

The discussion here is devided into two parts 
PART A consists of 

(1) A perspective view of water in clay- 
water system and its role in the 
abnormal flow throi:igh clays 

(2) Models describing the nature of bonded 
■ water 

PART B consists of a discussion about 

(1) viscosity 

(2) density 

and (3) thickness of bound water v/ith special 
reference to its role on the observed non-Darcy 
behaviour. 

PART A 

(1) A perspective view of water in the clay -water 

system and its role in the abnormal flow through 
clays, 

Karadi and Fagy (1961) and Williamson (1951) 
gave an excellent picture of the state of water in the 
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clay-water system. Water consists of polar molecules 
formed by a negatively charged oxygen ions and 
positively charged hydrogen ions. As the dielectric-,, 
constants of water and of the mineral, particle differ 
appreciably from each other, an electric field having 
surplus energy develops on the surfsice of the mineral 
particle, and as a consequence, water molecules become 
polarized and dipoles of the water become oriented. 
These oriented dipoles stick to the mineral particle 
surface. The range of electromole cular forces is 
around 0.25 to 0.5 and their magnitude sharply 
decreases with increasing distance from the surface of 
the soil particle. The magnitude of this force on 
the soil’ particle surface may attain a value of 10,000 
kg/sq cm, but at a distance of 0.5/Uit is neglegible. 
Citovies (1951) schematically presented this situation 

aiyi p -...i n the fig. the fluid of 
thickness 0.1 to 0. 15 X<. around the soil p'-irtiiefte 
adheres to it with a large force and hence it is 
called firmely bonded or adsorbed water. Beyond this 
layer the electromolecular forces decrease -appreciably, 
and the water molecules are .loosely bonded, sometimes 

reffered to as lios orbed water. This loosely 

bound water may exist up to an appreciable distance 
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and in fact the range of void ratios found in soils 
for engineering use are such that there is a very good 
probability that allmost all the y/ater there, is 
bound (loosely or strongly). In support, Grim (1952) 
writes that in a concentrated clay-water system as in 
a plastic^day ir&ss„ the nature of the adsorbed water 
is of' extreme importance in determining its properties 
because all or almost all of the water is non liquid. 
Macey (1940, 1942) also supports Grim’s idea. lutz 
and Kemper (1959) obsei*ved the same feature with 
Ka-clays. From vapour pressxnre data Shereshefsky (1928) 
as quoted by lutz and Kemper (1959) observed that the 
wall surface of a small capillary exerts an appreciable 
influence on the water molecules, even at a distance 
of 10,000 2 or about 3000 water molecule diameters. 

After a suxvej?" of the existing literature Romanov 
(1968)' summarized the pecularities of bound- water — 
as follows. 

(1) Bound water does not freeze even at very low 
temperatures , as low as - 78oc (Boioyancas 1917 and 
Rode 1952 ) 

(2) It’s density and -viscosity are different from 
those of free water, 
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(3) It’s specific heat is less than mity, about 

0. 76 cal/gm, i.e,, closely approaching that of 
ice (Andeanov 1949) 

(4) ’ lilms of adsorbed water possess shear elasticity 

1. e., Approach the mechanical 

properties of solid. 

(5) Although it is a non conductor, its dielectT*ic 
constant is smaller than that of free water 
(Dumanskii & Chapek 1934) and it is incapable of 
containing solutes (Dumanskii & Dumanskaya 1934). 

Martin (I960) also presented an excellent review on 
the different aspects of adsorbed water in clay. 

Macey (1940, 1942 and 1948) determined the change in 
permeability with v^ater content for saturated clays. He 
found that observed permeability values were much smaller 
than the calculated values from equation (5)- Macey 
observed that better agreement between the experimental 
and calculated values was obtainable by assuming that a 
very thick layer of bound ws,ter exists an the clay 
surface, 

Michaels and Lin (1954) conducted experiments on 
kaolinite bed and their results (fig. 7) show that the 
permeability decreases with increasing polarity of the 
permeant. The plots indicate that the Kozeny-Carman ■ 
equation is obeyed over a considerable void ratio range 
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for all the permeants. However, deviations from 
Kozeny -Carman relation occur for all permeants at low 
gradients and none of the curves pass through the origin. 

(1) Abnormal fluid properties (e.g., increased 
viscosity, changed density etc.) contribute less 
than 30 percent towards the total resistance to 
fliiid flow “through kaolinite beds. 

(2) Clay fabric changes during permeation with some 
permeants playing a predominant role in the non- 
Darcy behaviour at low gradients, 

(2) Models describing the nature of adsorbed water; 

Precise nature and structure of bound -water is 
not yet known, but it is generally agreed that some form 
of organization of the water molecules exists in the 
water immediately adsorbed on the clay mineral siurfaces. 
(Grim, 1962). 

At present there are three models available 
for adsorbed water structure and they are 

(1) "Ice like" or "Solid" structure proposed by 
Rosenqvist (1959) 

(2) Hexagonal net structure proposed by Hendricks 
and Jefferson (1938) 
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or (3) Two dimensional fluid model proposed by 
Martin (I960) 

Each of the theories is briefly described below, 

Eosenqvist*s "Ice like” or "Solid” structure 
hypothesis presumes that the water molecules are not 
rigidly fixed into^ a well ordered lattice but there is 
more "order" than in normal liquid water. One may infer 
that each bond is stronger than in normal liquid water 
due to the polarization of the clay surface. Presumably, 
this water would resist both normal and shear stresses 
to a larger degree than normal water. 

Hendricks and Jefferson’s Hexagonal net structure 
hypothesis states that the hydrates of the silicate 
minerals in montmorillonite , contain layers of water 
molecules joined into a hexagonal net. Purther they 
explain the stability of the hexagonal net by the attraction 
between hydrogen atoms of v/ater molecules and neighboring 
oxygen ions of the silicate layers, or oxygen atoms of 
other water molecules in the net. 

The two dimensional fluid model proposed by 
Martin (I960) recognizes a stronger bond in adsorbed water 
than normal water. Martin visualised the water to be 
more tightly bonded laterally than normally, resulting 
in the' ability to withstand large normal stresses but 
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not any shear stresses. Michaels (1961) supported 
Martin’s two dimensional fluid model. 

P AST - B 

('^) Viscosity of Water in Clay-Water System ;- 

A good amount of quantitative and qualitative 
data is available for the viscosity of loosely bound or 
strongly bcutd adsorbed water in clay water system. All 
the research T^/orkers agree that a water structure, which 
varies in extent with particle arrangement and the 
adsorbed cationic species, exists at the surface of clay 
particles r and this gives. a high viscosity to the 
adsorbed water, 

Rosenqvist (1961) observed that the viscosity of 
clay water increases with decreasing distance from the 
mineral surface. He indirectly determined the coefficient 
of viscosity from the coefficient of self diffusion 
data for different moisture contents and the results 
show an exponential increase of average viscosity with 
decreasing moisture content as reproduced in fig. 8. It 
is seen that an average viscosity of 24 centipoise is 
obtained at 30?^ water content and 153 centipoise at 10^ 
water content in contrast to the viscosity of 1 
centipoise of ordinary water not acted by any surface 



forces. Kemper, Maasland and Porter (1964) estimated 
■viscosity of water adjacent to sodium and calcium 
saturated bentonite surfaces from diffusion rates of 
deuterium hydroxide in oriented clay pastes at several’ 
moisture contents. Their data indicates aroimd 60 % 
reduction of mobility (inverse of viscosity) in water 
which is more than three molecular layers from clay solid 
surf ace . 

low (I960) calculated the viscosity of water in 
the clay water system from the activation energy data, 
using the relation 

,e/rt 


B 


(7) 


where » coefficient of viscosity 
0 B s= an empirical const. 

^ = exponential 
E '= activation energy 
E = Molar gas constant 
T = absolute temperature 

low estimated the activation energy E from permeability 
tests in conjunction v/ith equation 


log Q — log ( Ti ■^t) 


B 


1 _ 

ED 


>( 8 ) 


where K = Intrinsic permeability == 
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Q ~ llow rate 

A = Area of flow 

i = hydraulic gradient . 

Equation (8) is obtained from eqn. (7) and Darcy’s law 
(Q = K Ai). Equation (8) suggests that the plot of 

■j 

log Q vs ^ permits computation of activation energy (if 
the pressure gradient is kept constant), low’s data 
indicates that the activation energy for viscous flow in 
the clay-water system is approxanately 13?^ larger than 
the activation energy for viscous flow of pure water which 
clearly indicates a very high viscosity for water in the 
clay -water systems. Under similar conditions, however, 
Rosenqvist (1959) observed a 250 percent increase, like 
Rosenqvist (I960, 1961), Olsen (I960) also concluded that 
viscosity decreases exponentially with distance from the 
clay particle . One of the main reason why exist ing 
theories of hydraulic conductivity (like Kozery Carman’s 
equation) do not give correct results in case of surface 
active soils is perhaps the existing equations do not 
take this increased viscosity and its variation into 
account. Carman (1939) and Macay (1942) concluded that 
existing theories of v/ater flow through porous media can 
not be applied to clays or surface active soils without 
assuming that part of the water which is immobile to have 
high viscosity. 

V . 



27 


To study the effect of this increased •viscosity 
of strongly and loosely bound water on the flow rate 
or permeability, Olsen (I960, 1961) presented a capillary 
model consisting of infinitely viscous liquid of thickness 
’t* clinging around the pore wall, and beyond the distance 
’t’ from the wall the liquid is assumed to have constant 
viscosity equal to the bulk liquid value. This assumption 
simply means that the influence of high viscosity on flow 
rates is equivalent to a reduction in the pore size by 
a thickness 'f. Olsen (I960) says, that a mor^ reasonable 
concept of abnormal viscosity probably v/ould be one in 
which viscosity decreases exponentially with distance ' 
from a clay particle surface. Nevertheless, the use of 
simpler rigid film concept in this analysis can be . 
justified on the grounds that both concepts lead to 
similar violations of Poiseuille’s law. Using rigid film 
concept, Olsen derived the following flow rate equation 


^ozoni- 


= 1 + 


ksyx-^) 


,3 


1 


n" 1 

1-n^ T^ 


(S) 


where there is a rigid 

film of infinite viacDuo l.-yer around 


, clay p article 
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%02:en7 “ rj^to for Koaenjr Cr-ri'ian'B cqv.e.tion 

with constant viocoBit^'" of pore wo-tor 
= Specific surface per unit volune of porticles 
t = Thickness of rigid f iln 
n = Porosity 
T = Tortuosit 3 r factor 
K© = Pore shape factor ^ 2.5 

Equation (9) gives the estiir-utsd possible cLiscrepencieo due 
to high viscosity in water of clay -water* system, Eqn. (9) 
is plotted and is shown in fig. (3) • Per interpret at ion 
of this fig. the assrimption needed is that *t’ reo-xins 
constant in the entire flo\.' process. Equation (S) ^ 

fig* (5) shov? that for po.rticulor thicknoss ’t*^ qb porosity 
increases discrepencj^ ratio (= approaches 

unity but at lov/er porosities, nctual flow rote 
is nuch smaller than preSictoci hy hosenjr— Co.r:uc:,n’s equation 

^%oz.Qny^* 

. Blaciaaore (1963) experiEiantallx'- confirBied the 
effect of adsorbed v.ater fill', thic’raess s.nd its viscosity 
on the hydraulic cOirductivitxr of bentonite and kaolinits. 


(2) Density of Adsorbed or Bound 'fater:- 


Eealising the importance of role of bounci -water in 
the physical properties of clay water systen, several 
investigators tried to find the density of bound water. * 



CaJfimXated Atitaiuj j from Stirfac« in 
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Although there is qusutite-tive difference between 


the VLilues of density of adsorbed v/ater found by various 
sutiiors but qualitatively?- it is observed by IJorrisb. (1954), 
Anderson & lov/ (1958), Girahar:: (1S?62) a,nd Dohney and 
IT'ungairoli (1972) that density of water in clay?" -water systei^is 


above 5,7; water content is less than 1,0 gn’/cc (nii-niiaun 
vaiue quoted is 0.97 gr/cc) and as the viratsr content 


i-ncreases beyond that the density 
1,0 gir/cc (density of free water) 


gradually e.pproaches to 
ihepariaionts nerforued 


by :De’.;it md Arens (1950), llachensie (1958), Kooney et al 
(1951) indicate that "oqIoxj 5;-5 vi'ater content the density 
shaxplyr increases beyond 1.00 gro/cc end reaches a peak 
value of 1,4 gin/cc indicating;; thereby an increase of 
density as one approaches the clay surface, The results 
of Anderson & low (1958) are shown in fig. 10 and table 1. 
Density data for Dewit and Arens (1950), Forrish (1954) and 
Mooney et el (1951) are shown in tables 2 & 3. Dor better 
perspective view about the density data covering 

a wider water content range, Dohney & Dungaroli (1972) 
plotted the results of all research workers including 
themselves and is shown in fig,. 11 * 

Buchtnoom (1964) observed an increase in the densi-tgr 
of water (from an average of 0.98 to an average of 1,13 
gir/cc) in a consolidating soil at the end of consolidation. 
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The increase of density with. the reduction of water (or 
closeness of the soil particle) is attributed to the fact 
that more and more saturating water falls within the surface 
forces of the solid minerals and thereby/ making the water 
more and more abnormal compared to the free water, 

Kutilek .(1964a, 1962) observed that the density 
of adsorbed water is lower than 1-oand is influenced by the 
surface properties of the minerals and by the exchangeable 
cations. He further states that water density is higher on 
the MontoMJorillonites than kaolinites because monovalent 
cations cause higher density than bivalent cations. 

(3) Thickness of Bound Yfater; 

The thickness of adsorbed or bound water is of . 
vital importance in the quantification of total probable 
discharge from surface active soils under hydraulic gradients. 
According to G-rim (1953) many attempts have been made to 
estimate the thickness of adsorbed water v/ith definite non 
liquid characteristics without very good agreement among 
the values suggested. HouwihLc (1937) arrived at a value of 
25 . G-rim & Cuthbert (1945) suggested the thickness of 
strongly adsorbed water to be 7.5A° to lOA® with'gradual' 
transition to liquid water. Jaeger (1938) and Spiel (1940) 
concluded that the thickness of adsorbed water at optimum 
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TABIE 1 


DEISITY OP ADSORBED 

WATER AT 25 


(Data from Anderson & Low (1958) 

(1958b)) Quoted in Martin (I960) 

and Low and 

Anderson 

Water content 

Adsorbed water density (gir/ml) 

(gm HgO/gm clay) 


Li 

K 

0.60 

- 

- 

0.980 

0.78 

0.9715 - 

- 

- 

1.00 

0.9731 

- 

0.983 

1.214 

- 

0.977 

- 

2.00 

0.980 

0.984 

0.988 

3.00 

0.986 

0.989 

0.993 

3.35 

- 

- 

0.994 

4.00 

0.991 

0.992 

- 

4.61 

- 

0.993 

_ 

4.80 

0.995 

- 


6.70 

0.9971 

0.9971 . 

0.9971 
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• TABLE 2 

ABSOHBEB WATER DENSITY OK CLAY 


( Dewit & 

Arens (1950), 

Quoted 

in 

Martin (I960)) 

Mineral 

Relative 

Humidity 

(percent) 

Water content Density of 

Adsorbed water 

(mg/g) 

Montmor illonit e 

25 

116 


1.41 


50 

165 


1.37 


75 

284 


1.32 


100 

460 


1.16 


- 

2440 


1 .02 


- 

3010 


1.02 

Kaolinite 

25 

4 


1.68 


50 

8 


1.12 


75 

9 


1.03 


100 

65. 


0.99 

Illite 

25 

30 


1.36 


50 

49 


1.17 


75.' 

69 


1.08 


100 

188 


1.04 



1038 


1 .00 
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TABLE 3 

BEBSITY OB ABSOEBEB WATER GALGUIATEI) EEOM X’RAY DATA 


(Data from Mooney, Keenan & Wood (1951) and Korrish (1959) 
Quoted in Martin (I960)) 


Author 


Experimental Data 


^100 

(A) 


Y/ater content 
(gH 20 /g clay) 


Calc\ilated adsorbed 
water density 
(gm/ml) 


9.8 

12.4 

12.4 

Money et al 15.4 

15.4 

9.5 
19.0 
40,0 
61.0 

73.5 


0 


0. 10 

- 

0.15 

1.46 

0.20 

- 

0.28 

1.27 

0 

— 

0.37 

0.986 

1.19 

■ 0.986 

2.0 

0.986 

2.5 

0.99 


Korrish 





plasticity of kaolinite is of the order of several hiondred 
angstroms. Rutgers (1954) suggested a value of 16 to 
lOOA® . lutzor & Kemper (1959) calculated the thickness 
of water adsorbed on Wyoming Bentonite siJirface from the 
p2?essure memlrave data for different equillibrium pressures 
and the results are shown in table 4. Norton & Johns-en-" • " 
(1944) reported the thickness associated with kaolinite 
particles ranged from 2901° to 13A° with increase of 
pressure, Williamson (1951) reported a thickness as high 
as 1000A° which means allmost all the water in the natural 
soils within the working void ratio range is bound water 
(loosely or strongly). This is also supported by Grim 
(1952). 
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TABLE 4 

Thickness of- water layer on surfaces of clay particles 
(Wyoming Bentonite) 

(Ajfter Lutz & Kemper (1959)) 


Thickness* of water films at 
5 equillibrium pressures 




O 

A 

<5 

A 

o 

■ A 

1 

O' 

A 

Ka-Saturated clay in water 

66.0 

35.9 

25.0 

15.1 

9.4 

Ka-Saturated clay in 

0.03 K laCl 

61.5 

33.3 

24.2 

15.5 

8,5 

Ka-Saturated clay in 

1.0 K NaCl ■ . 

46.2 

24.1 

16.0 

10.3 

6.8 

Ca-Saturated clay in water 

40.5 

20.^ 

14.8 

9.6 

7.1 

Ga-Saturated clay in 

0.01 K. CaClg 

35.8 

18.9 

14.4 

10.6 

7.3 

K-Saturated clay in water 

68.2 

36.6 

24.5 

14.3 

8.6 


♦ Calculated on the basis of 800 m^ of surface/gm, of clay, 

' 

and 10 A as the packing area for a single water molecule. 
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Chapter 3 

EXPERBlEITjtt MID THEOEETICAE FOR DAECY VELOCITY 
GEALIERT EELATIOISHIPS 

Various authors have suggested several empirical 
velocity gradient relationships based on their experimental 
results. There are also a few theoretical solutions available 
which are predominantly based on capillary tube model.. 

Kutilek ( 1969 ) suggested several v-i relationships based on 
non Rewtonian flow. The types and names of equations will 
be given first and their discussion will follow. The 
chapter is devided into 4 Parts as given ’below 
'A*: E!£PIEICAL IQTJATIORS AVAILABLE 
•B': THEORETICAL SOLUTIORS AVAILABLE 
*C's SOLUTIONS BASED ON NON NE'^,T?ONIAN ELUID FLOW 
•D’ ; SHAPES OE EXPERIMENTALLY OBSERVED VELOCITY^ 
GRADIENT CURVES. 

PART A . 

EMPIRICAL EQUATIONS AVAILABLE 

The following is a list of empirical equations 
reported in literature. 

(1) Puzyrevskaya’ s equation (1931) 

V = K(i - i„) 


( 10 ) 
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(2) Valarovich and Tchuraev's Equation (1964) 


■V = Ai + ^ - B 
i 


(3) Izbashe's Equation (1931) 

.n 


( 11 ) 


V = Ki“ (12) 

(4) Slepicka’s Equation (1961) 


n f f 

V = .o<(|) 1 


(13) 


(5) Swart zendruber’ s Equation (1962a5 1962b) 

I”- 

i) V = M i - I( 1 - 4 ; 

ii) V = B 1^1 - J(1 - 

«Q<d 

(6) Kutilek's Equation (1965) 

V = M 


— 

..... (14) 


(15) 


i log (A + e®^) - 1^ 


(16) 


(1) Puzyrevskaya ' s equation ; 

Puzerevskaya (1931) first proposed the equation of 
the form 


V = K(i - i^) 


10 


where i^ =■ Threshold gradient below which 


there is no flow. 
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Li ( 1963 ) found that for Houston and varved clay. Threshold 
gradient decreases approximately linearly with increasing void 
ratio and his results are shown in fig. 12. Kondon (1967) also 
observed that the threshold gradient increases with decreasing 
porosity and increasing specific surface and he explained his 
results on the basis of rheological properties of water in thin 
slits and narrow capillaries. Bondarenko (1968) concluded' that 
liquids with H-bonds behave as Bingham material and have a 
definite threshold gradient due to certain yield stress of the 
pore water. About the validity of equation (10) "Valarovich 
and Tchuraev (1964) suggest that the equation gives correct 
result only when one works with large pressure gradients and 
materials with uniform pore radii. 

(2) Yalarovich and Tchuraev* s Equation. 

Based on experimental results, Yalarovich and 
Tchuraev (1964) forwarded an equation of the form 

V = Ai + ^- B (11) 

where A, B and C are constants which depend on 

the pore radii, the type of pore distribution, 
the porosity of the sample and the rheolo- 
gical properties of the fluid. 

This equation is identical with Herpin's' Bingham body equation 
to be discussed later. 
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(3) Izbash's Equation; 

Eor the non linear portion in the laminar region, 
Izbash ( 1931 ) proposed the equation of the form 

V = Ki^ (12) 

where n ^ 1 depending on soil properties. 

Hans bo (I960) supported the same equation and claimed 
that his data and Silfverberg' s (1949) data fit well in the 
above equation for the values of 'n* in between 1 to 1,5. 
Dudgeon (1966) found the value of n to be 1.12 to 1.25. 

(4) Slepicka's Equation; 

Erom dimensional analysis and experimental results, 
Slepicka ( 1961 ) proposed an equation of the form 

V = CX(|j kV (13) 

or V = i^ (13a) 

. f-1 ^ 

where = c<(^ K 

= general coefficient of filtration 
^^ = Dynamic coefficient of viscosity 
q— = Surface tension (contact tension of 
solid liquid interface.) 
cx^= A factor. . 
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Depending on the values of f in equation (13) or (13a), the 
flow regime can be devided into preline ar regime, linear 
regime and post linear regime as explained in fig. (13). 

Slepicka observed an experimental linear relation 
between the value of f in the preline ar regime (our main 
concern) and Darcy's K (slope of v-i curve in the linear regime) 
in log scale as shown in fig. 14. Slepicka also found a linear 
plot, between critical gradient i^ (as defined in fig. 13) 
and Darcy's K. The result is shown in fig. 15. As seen from 
fig. 14, the value of f (which is nothing but 'n* in Izbash's 
eqn. Eq. 12) varies between 1.3 to 2.3 depending on Darcy's 
K value. 

(5) Swart zendruber ' s Equation 

Erom the experimental v-i curves Swart zendruber 
(1962a) proposed the equation 

V = M 

where M and I are constants. I is the' 'i' 
intercept obtained when linear part of equation 
(14) is extended as explained in fig. 16. 

I is the measure of departure from Darcy's law. 

The fit' of the equation (14) for the experimental 
data of Hansbo (I960), Von Bngelhardt and Tunn (1955) and 
Lutz «Sb Kemper (1959) are shown in fig.- 17, 18 ^6, 


1(1 - 




(H) 





Linear regime 
f=1.0 


Post linear 
regime ,f^lo 


,k 


.Pi’o^ine r Re gime , .0 


irlyc: raulic Gradient 

FiF,.,i;j; Velocity Gradient Relatlori-,;hli. as pr'-.M:’v^s...a 


.cka ,1961 




Fig. 14s Hci lati on between 

f k (Pnj-linour) 
( » ’ ' I; 'j r . j 1 e ! c v: : * , j > ] ) 


lO-^ 10’ 10' 10 

K » 

Fig. 15 » Relationship between 
liA B: 

(After Slepicka,i96l) 





Superficial Velocity iv) 
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Hydraulic Gradient 


Fig. 16: Velocity Gradient Relationship as proposed 
by Swart zend ruber (1962 a) 






Flow '/elocltyjlO'-^Cin/Se 


‘-f i 


I 

1 . 

■ 



Hydraulic Gradient 

Fig. 20 s Plots of equation 15 with B-“i & J - i 



Hydraulic Gradient , Water Basis 

Fig.21i Bita for Na- sat. Utah Bentonite (lAitz & Kemperfl959) as fitted! 
. by eqn.lS (Solid line ) and eqn. 14' (Broken lines) I 

After Swart 26 nd ruber, 1962 b I 
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Later, Swart zendirLter (1962b) improved equation (14) 
by bringing one more parameter and his equation reads as 



J(1 


€ 



(15) 


where B, J and C are constants. 

Equation (15) is plotted in fig. 20. The term J^C is a 
measure of non Darcian behaviour. 

Swart zendruber claims that equation (15) fits better 
with the experimental results than equation (14) mainly because 
in Bq. 14 tho slope of the v-i curve at i = 0, is always zero, 
regardless of the magnitude of I, (except when 1=0). A 
'comparison of equation (14) and (15) for Lutz and Kemper’s 
(1959) data is shown in fig. 21. 

(6) Kutilek's Equation 


Assuming a higher viscosity of pore water near clay 
surface and based on experimental results Kutilek (1964) 
proposed the equation 

■"1 

( 16 ) 


V = M 


1 log (A + e®’-) - I 


0 


♦ » • • » 


M 

los mT 


whe re 


B = 


A 


M’ 


- 1 


M = hydraulic conductance 



m 



'vHl 



A*jCnaot0A 


Fig. S3' Velocity Gradient Ifelatlonsh 
ais proposed fc-y i'srpin i^nd 
Tchudnov ski j <13675 
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M' = Initial hydraulic conductance 
= Intercept of M on I axis 

The parameters M, M’ and I^ of equation (16) are explained 
in fig . 22 . 

Typical values of M, M' and I^ for different soil 
mineral -liquid systems are shown in table 5. 

PAST ‘S' 

THEORETICAL SOLUTIONS 

The following theoretical solutions available 

(1) llerpin's Bingham body equation (1967) 

(2) Nerpin and Deryagin's equation (1969) 

(3) Kovack's equation (1967) 

(4) Klausnor and Kraft's equation (1966) 


(1) Nerpin's Bingham body Equation; 

Nerpin and Tchudnov ski j (1967) postulated that soil, 
water behaves like a Bingham body and siiggested an equation 
identical with the flow of Bingham liquid through a system 
of capillary tubes which reads as 


V = 


K I 
o 


3 
U 




f + 1 


( 17 ) 
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where 


/ S 

T = -2 — = threshold gradient 

0 n V 

T = Threshold yield stress 
o 

3 = Specific surface 

. n = porosity 


Y = Density of water 


Eq = 


8 


hydraulic conductivity for a 
Newtonian liquid 


r = radius of capillary 


y= Density of liquid 

Viscosity of liquid. 

Equation (17) j takes the shape as shown in fig. 23. 

Bondarenko (1968) found the value of yield stress (T^) 
of the order of 10“^ to 107 ^ dynes/cm^ for liquids with 
hydrogen bonds e.g. H 2 O 3 C 2 H^ 0 H, CH 2 OH etc, 

(2) Nerpin and Deriyagin's Equation; 


Nerpin and Deriyagin (1969) derived a flow equation 
from thermodynamic principle as 


9 = ^ grad 0 + Ss^ady ... . ( 18 ) 

where 0 = moisture potential 
>^ ■= viscosity 

i.i.T. KANFUR 
^FNTRAL I JBRAgY 

Tsias 


Act., 
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TAEEE 5 

pARil'-IETERS M, M' and FOR DIFF . CMY MIHERAIS 
(After Kutilek 1965) 


Cation Exchangable 
Mineral 

t 

liquid 1 

? 

! 

! 

f 

L. 


MxlO"^ 

cm/hr 

! 

; M'xio"^ 

1 cm/hr 

t 

f 

t 

H-Kaoiinite 

H^O 

160 

23.7 

1 .5 


OtOI M HCl 

35 

30.0 

19.3 


0.001 M HCl 

105 

24.6 

8,3 

Na Eaolinife 


395 

8.1 



Ofl H HaCl 

0 

22.2 

22.2 


0.01 n Had 

0 

19.5 

19.5 


0.01 H Ha2S0^ 

202 

15.5 

3.3 

Ca Kaolinite 

dO 

11 5' 

14.8 

1.5 


Oil H CaCd 

0 

24.9 

2f.9 


0.01 H CaCIp 

0 

18.6 

18.6 


0.1 I Ca/TOA 

20 

16.4 

0.2 


0.01 H Ca/HO^/g 

113 

14.8 

1 .5 

H-Montmor illoni t e 

HpO 

171 

1 .60 



OfOl H HCl 

110 

1.93 

0.85 


0.001 I'l HCl 

138 

1.45 

0.37 

Na Montmorillonite 

HpO 

340 

■ 0.'35' 

UTUj 


1 H Had 

97 

• 31.2 

5.3 


0.1 H Had 

255 

12.7 

1.13 


0.01 Had 

325 

8-4 

0.10 

Ca Montmorillonite 

HpO 

295 

4.20 

0.12 


1 H CaClo 

90 

111 

20 


0.1 H CaCd 

120 

108 

9.3 


0.01 H CaCIj 

210, 

66 

4.3 

H-lllite 

— 


18.2 

1.1 


ofl H HCl 

140 

31.4 

4.2 


0.01 H HCl 

160 

24.5 

1.5 


0.001 H HOI 

260 

88 

1.1 

Ma-Illite 

HpO 

1B5 

" OTBT 



1 H HaCl 

52 

3.6 

0.64 


0,1 H Had 

85 

1.21 

0.38 


0.01 H Had 

110 

0.93 

0.33 


0.1 H HagSO^ 

100 

1.02 

0.38 

Ca- ILlite ./ 

__ 

250 

13.2 

1.5 


1% CaCd 

150 

22.5 

4.2 


0.1 H CaClo 

185 

15.4 

2.2 


0.01 H CaClg 

210 

13.8 

1.5 ; 


1 
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S = Specific kinetic energy 
^ = Active porosity 
0(' = A nnmericel coefficient 
S = Coefficient of convolution of the 
pores 

This is ary one of the following three 
quantities 

(i) Temperature 

(ii) Electric field potential 

(iii) Concentration of dissolved 
substances 

Corresponding coefficient of thermo-, 
osmosis or capillary osmosis 
After making the following substitutions 

■3 . 


K 


?( £' 

grad 0 = i 

.2/ . r2 


= 1 , 


- MS grady- 

equation (18) gets converted to 


q = K(i - Iq) 


• . • • « • 


(18a) 


(3) KOVAC's EQUATION. 


Based on capillary tube model and with liquid obeying 


the relation 
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r= r 


L Jn 


(^9) 


where 7 „ = Threshold shear stress 

. 0 

= Viscosity of fluid^'n' symbolises the 
direction^ at right angles to the 
direction of movement 

Kovac (1967) derived the following velocity -gradient relationship 

. 2 


V 

T 


D 


(1 - ^) 


2 to 

3 i 


ir, 3/2 

(1 - tan / 1 i 


log/-(1 -/) 


. . .. . (20) 


where v^ = Darcy velocity for Newtonian fluid 


8 




Drom his experimental results, Kovacs approximated equation 
(20) as 

i - i 


V = 0.714 ( 1 - -) 


( 21 ) 


for i <^12 i^ 


where = Darcy's permeability coefficient 




r 


Eqn. (20) and (21) axe plotted in fig. 24 and 25 ^ 


(4) Klausner and Kraft ' s Equation; 




Based on capillary tube model Klausner and Kraft 
(1965a, 1965b, 1966a and 1966b) theoretically explained the 


Plg»^t ai3*tt<mshlp l»tween 7/Ei^nd 1/i 
(After KCTr«es, 1967) 
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non. Darcy behaviour of soil water flow assuming wall forces 
arising out of active clay surfaces. They considered a force 
P acting normal to the wall as shown in fig. 26. The force 
P is the average electrical field arising from the unbalanced 
charge on the crystalline wall. This force ’P' is the same 
electrical force that gives rise to Cihapman-G-ouy potential 
at a crystalline interface. The force P is assumed to influence 
the fluid by bringing into existence a force f of magnitude 
, C is a const, of proportionality) acting along the 
capillary axis (see fig. 27) and resisting the motion. 

Fort a, bvi eo^t cU / <5^ -n-ecx/t -tk-x. 

and Kraft 

had set up equillibrium momentm equation for a Newtonian 
fluid and then solved for discharge Q for different pressure 
gradient (/ii J ) and the ultimate solution reads as 

4 r~ S' 

- Q(Aj, R) = |- AJ - -I a^ W(a, R) - r^f(r, R) dr 

1 — o 

J 

Aj - R*(Ad, R) . , . .(22) 

where K|^ = Viscosj.ty of the fluid 

a = a(Aj 5 E) and is defined by the radius 
'r' upto v/hich f low zone extends and 
this is a solution of equation 

a(^) = W(a, E) where 
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I 


1 

2 

Direction of Flow 

g «'?■'*! W'TII force as assi^mej 
b- Kiausnyr and Kraft 
C19G6 b) 


2 



Fig*27* Force opposing the 
fluid motion hs 

V 1 11 ''C d I'V K 

brail ■ 
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r 

W(r, R) = ^ ^ '£{^9 R) dr 

o 

Eqn. (2-2) is shown in fig. 28. 

Equation (22) as seen from the plots qualitatively 
explains the observed non Darcy velocity -gradient relation- 
ship - but one of ' the main drawbacks of the theory is that, 
it does not consider the variable viscosity, density and non 
Newtonian behaviour of the pore fluid whose existence is more 
or less certain due to the very nature of clay-water 
environment. So an approach similar to Klausner and Kraft 
considering the soil-water interaction and rheological 
parameters of pore fluids holds a promising future. 

PART C 

SOLUTIONS BASED ON NON-NEWTONIAN FLUID FLOW 

Kutilek (1969), Govier and Aziz (1972) and others have 
suggested a few velocity-gradient (v-i) relationships based on 
non Newtonian fluid flow. Each non Newtonian fluid is charac- 
terised by a distinct shear stress (T) - velocity gradient 
(~) relationship. Ten v-i relationships can be suggested 
based on ten different non Newtonian fluids characterised by 
respective constitutive equations as completely described in 
Table 5 A. This table also gives velocity distribution and 
discharge through a capillary tube for some of the cases. If 
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uhapu of Discharge Pres^iure Gradient curve as 
found by Klausner and '.raft (l966b) 
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desired, for soils which can be modelled as a bunch, of 
parallel capillary tubes, the tube radius E can be equated to 
the soil porosity E and its specific surface S by the 
relation R = — g (Kutilek 1969). Eor the sake of general 
interest, a typical derivation of a v-i relationship for 
pseudo -plastic fluid is, given here, v-i relationship for 
other fluids can be derived similarly. Table 5A gives only 
the final result. 

A TYPICAL DERIVATION; 


Pseudo Plastic Pluid ; 

The constitutive eqn. is 

j n 

Z = m S = m( - ^) 


T ut d."? d*P 

= “ 2 dp ’ dL press, drop in a length 

of dL. 

Combining the above two equations, we have 

^ _ ( 1 - 
dr “ ''2m dL^ 


Integrating this equation and substituting u = 0 at r = R, 
we have 


u 



J_ dP 
2m dL 


1/n 


) 


n 

1+n 


R 


1 4~n 
n 


t+n 

n 
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P/iET »D' 

EXPERIMENT ALLY OBSERVED TYPES OF NORDAROY VELOCITY GRADIENT 
CURVES 

Kutilek (1969) and Olsen (1965) made a survey of exper- 
imentally observed velocity -gradient curves by various authors. 
The net results are summarised in figs . from (29 to 35). Each 
fig, gives the shape of v-i curves observed along with the 
names of authors, type of soil and moisture condition. 
Deviations met by various authors at low gradients are only 


considered here. 
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Chapter 4 

REVIEW OE EXPERIMERTia DATA 

This review is done under the following heads. 

A. Eon Darcy flow data for clays and clay minerals 

B. Eon Darcy flow data for sands and other inert 
materials 

0. Eon Darcy results attributed to causes other than 
the surface phenomena 

D. Summary of published data and analysis of 
deviations from Darcy behaviour. 

Results are reported chronologically as far as possible. 

A. EOE DmOY BLOW DATA EOR CIAIS AED CLAY MIEERALS 

One of the earliest non Darcy flow rate -gradient data 
for clay is due to Stearns (-1927) as quoted by Swart zendruber 
(1962b)* Stearns provided data on clay content ( <('5M ) 
along with the velocity gradient data and Swartzendruber {1962b) 
fitted equation la to the above data 

q=C^i+C 2 .(!a) 

where q = Blow rate 

i = Hydraulic gradient 
C^ and Cg are constants. 
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Fig.36» Gradient intercept i * Vjg* Content of tising 

Steam’s (1927) data (After Swart zend ruber, 196 2b) 



Jljd i'au.li 0 


Fig. 37* Results vif perffleability tests on muddy clay from 7asby 
perforiWd by Silfverb^' rg in 1949. Sample consolidated 
under a pressure of j.3 Kg/Cnr Curves drawn by Hansbo,l960 
(Aftf r Hansb''^}l960) 
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by the method of least squares. Prom equation (la), the 
gradient intercept (i') for q = 0 is obtained by the equation 

i’ = ^ (2a) 

and i* can be considered as a measure of deviations from 
Darcy's law. Swart zendoruber calculated the i' values from 
Stearns' data and plotted them against percentage of clay 
( by weight as shown in fig. 36. 

Por the samples of low clay contents, the calculated 
points are bounded by i' = + 0.06 signifying thereby the 
obeyance of Darcy's law within reasonable experimental error. 
The remaining i' values outside the selected error interval of 
+ ©.06 are positive, and have clay contents above 4, percent 
Swart zendruber concludes that fig. 36 gives a strong 
indication that deviations from Darcy's law are associated with 
clay content which is a crude measure of surface active 
material. 

Silfverberg (1949) as quoted by Hansbo (I960) conducted 
some velocity-gradient tests on "muddy clay" at low gradients 
and his experimental results show a nonlinear relationship. 
However, no conclusions can be m§,de regarding the existence 
of threshold gradient as shown in fig. 37. 

lutz and Kemper (1959) studied the effect of water 
film (bound) on the intrinsic permeability of different clay 
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layer systems. The materials tested were Ma, Ca and H 
saturated Bentonite, Halloysite, \¥oyming Bentonite and Bladen 
clay. The percolating solutions used were? for Ua-clays, 

I'aCl; for K clays, KCl^l for Ca alays , CaOlg and for H clays, 
HOI. The normalities of the solutions were, varied . The 
experimental results are shown in Tadoles 6 , 6 A and 7. The 
experimental results indicate the increase in intrinsic perme- 
ability with increasing pressure which is completely out of 
accord with Darcy's law. This was p 8 .rticularly true for Ha- 
clays, in which the diffuse l^er of cations and the water ' 
structure extended to a distance greater than the radius of the 
pores between particles 5 this was attributed to a breakdown 
of the water structure. In Ca-clays with a restricted diffuse 
layer and water structure, increasing the pressure head 
particularly had no effect on permeability. 

Prom experimental investigation on Sweedish clays, 

Hans bo (I960) confirmed the nonlinearity in the velocity 
gradient relationship at low gradients. However, no threshold 
gradient was found and this nonlinearity was attributed to the- 
interaction between skeleten grains and pore water. The 
results are shown in Pigs, 38, §9 and 40. 

Karadi and Hagy (1961) hypothesised that the lower 
validity limit of Darcy law and the v-i relationship in the 
low gradient zone is controlled by rheological properties of 
bounded water sar rounding the soil particles. They further 
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TABia 6 & 6A 


(After Lutz and Kemper 1959 ) 


Effects of pressure, Svdsorbed cations and electrolyte concentration on 
the intrinsic permeability of clays. 


System 

P 

cms , 

of Hg 

Percolating Solutions 

- - . ' 

Water 

. _. ... 

i Sal-t 

- . _ 

or Acid* 

i Salt 

_,! - 

or Acid* 

-^xIO"^ 

cc/sec 

L, aq/ t 
r "Ap 

t 2 -S 

; (zio ^ 

) |cc/sec 

f 

! 

. . 4 - - - . , 

A P 

^(xIO"^) 

| A ^ 10 “^ 

{^it 

|cc/sec 

» 

1 

t. 

4. ^ sZ t 

! AP 

i ^(xio 5) 

! 

f 

f 

r 






( 0 . 005 N ) 


( 0 b 5 F ) 


10 

0.037 

3.18 

0.092 

7.81 

0.41 

35.30 

Efa-Sat . 

20 

0.089 

3.82 

0.230 

9.90 

0.90 

38.78 


30 

0.165 

4.76 

0.395 

11.30 

1 .43 

41.00 

Utah 

40 - 

0.241 

5.18 

0.610 

13 .10 

1 .99 

42.70 

bentonite 

50 

0.335 

5.75 

0.850 

14.63 

2.60 

44.80 


60 

0.460 

6 .60 

1.140 

16.36 

3.27 

46.90 


70 

0.610 

7.50 

1.470 

18.05 

3.95 

48.50 






(“ 0 . 005 U } 


( 0 . 5 U ) 


10 

0.66 

87.0 

1.32 

166.0 

2.50 

329 


20 

1.32 

87.0 

2.45 

153.8 

4.50 

296 

C a"“Sat • 

30 

1.90 

83.3 

3.58 

150.0 

6.50 

286 

Utah 

40 

2.50 

82.0 

4.65 

146.3 

8.50 

280 

bentonite 

50 


82.0 

5.65 

141.7 

10.5 

276 


60 

3.67 

80.5 

6.63 

138.6 

12.3 

270 


70 

4.33 

87.5 

7.55 

135.4 

14.0 

264 






( 0 , 005 U ) 




10 

0.72 

92.2 

0.82 

100.8 




20 

1.44 

92.2 

1.67 

102 . 2 



H-Sat . 

30 

2.16 

92.2 

2.41 

08 . o 



Utah 

40 

2.93 

93.8 

3.22 

98.8 



bent onit e 

50 

3.72 

95.2 

4.05 

99.3 




60 

4.57 

97.5 

4.93 

101.0 




70 

5.40 , 

98.6 

5.84 

102.0 








( 0 . 005 U ) 


( 0 . 5 U ) 


10 

0.68 

83.5 

0.90 

1 1 3 .0 

2.45 

322.2 


20 

1.45 

89.0 

2.0 

125.6 

5.0 

329.0 

Ua-Sat. 

30 

2.25 

92.0 

3.10 

129.7 

7.45 

W26.1 

halloy 

■40 

3.14 

96.5 

4.20 

131.5 

to . 10 

332.5 

site 

50 

4.05 

99.5 

5.28 

132 .6 

12.94 

340.0 


60 

5.05 

103.3 

6.32 

132.2 

15.79 

346.0 


70 

6.05 

106.1 

7.50 

134.6 

19.10 

359.0 


Continued . . . 
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Table 6 & 6A continued. 



10 

1.40 

168.7 

1.54 

“■"(■O'TDDWT 

193.4 


20 

2.85 

171 .7 

3.10 

194.5 

C 0*""S B'ij • 

30 

4.30 

172.7 

4.64 

194.0 

halloy- 

40 

5.75 

173.3 

6.25 

196.3 

site 

50 

7.20 

173.6 

7.80 

196.0 


60 

8.65 

173.7 

9.30 

194.7 


70 

10.10 

174.0 

10.90 

195.7 



10 

0.60 

86 .6 

0.78 

(0.0051T) 

98.0 

1.97 

(0.051T) 

125.0 


20 

1.44 ■ 

90.4 

1.58 

99.2 

1.96 

126.3 

H-sat . 

30 

2.24 

93.7 

2.34 

98.0 

2.94 

126.3 

halloy- 

40 

3.04 

94.1 

3.07 

95.1 

3.90 

125.8 

site 

50 

3.75 

94.4 

3.92 

98.4 

4.88 

126.0 


60 

4.55 

95.4 

4.73 

99.0 

5.86 

126.0- 


70 

■ 5.42 

97.1 

5.64 

101.0 

6.84 

126.0 



73 


TiSIE 6A 


i 

r 

i ^ 

SYSTEM ;cms 

; Hg 

1 

f 

f 

f 

I 

- - -1 - 

Percolating solution 

- ■ J 

Water 

Salt 



or Acid* ; 

Salt 

or Acid* 

Mx10~^!k' ^ 

^ 1 AP 

cc/sec 1 ^(xio”^) 

1 0”^ |k ' = — 1 

/it ; A p ; 

f , to ^ f 

cc/sec 1 (xIO ) ] 

f f 

^10-5 Ik 

cc/sec } 

T 

, aq/ t 

A p 
^(xIO"^.) 





(0.005IVT) 


(0.5F) 

10 

0.56 

0.69 

1 .13 

1 .39 

4.2 

5.90 

20 

1.26 

0.77 

2.48 

1.52 

9.6 

6.75 

Na-Sat. 30 

2.09 

0.85 

3.90 

1.60 

15.6 

7.32 

Bladen 40 

'3.01 

0.93 

• 5.35 

1 .64 

21.3 

7.50 

clay 50 

4.13 

1.01 

6.90 

1,69 

27.2 

7.66 

60 

5.35 

1.09 

8.50 

1.74 

33 .7 

7.90 

70 

6.60 

1.16 

9.90 

1.74 

39.7 

7.97 





(U.0U5K) 


(0.05KfT ~ 

10 

10.1 

13.03 

10.4 

12.76 

17.7 

22.21 

20 

20.0 

12.96 

20.8 

12.76 

35.6 

22.35 

Ga-Sat 30 

29.8 

12.80 

31.2 

12.76 

53.5 

22.40 

Bladen 40 

39.7 

12.63 

41.7 

12.78 

71.5 

22.45 

clay 50 

48.9 

12.60 

52.0 

12.76 

89.0 

22.37 

60 

■58.1 

12.50 

62.3 

12.76 

106.2 

22.25 

70 

66.9 

12.30 

72.6 

12.73 

123.8 

22.20 







(0.05if) 

10 

9.0 

10.60 

11.5 

13.88 

12.2 


20 

18.0 

10.60 

21.8 

13.14 

23.4 

14.38 

H-Sat. 30 

26.8 

10.51 

31.2 

12.53 

33.5 

13.62 

Bladen 40 

35.5 

10.45 

40.5 

12,20 

44.4 

13.65 

clay 50 

44.0 

10.27 

49.6 

11.96 

53.4 

13.13 

• 60 

52.4 

10.27 

59.0 

11.87 

63.6 

13.03 

70 

60.5 

10.17 

68.5 

11.79 

73.6 

12,81 





(0.005Jf.) 



10 

0.27 

0.31 

0.34 

0.40 



20 

0.51 

0.31 

0.67 

0.39 



Ca-Sat. 30 

0.76 

0.31 

1.01 

0.40 



Wyoming 40 

1.05 

0.32 

1 . 34 

0.40 



hento- 50 

1.32 

0.32 

1 .69 

0.40 



nite 60 

1 . 60 

0.32 

2.08 

0.41 



70 

1.89 

0.33 

2.53 

0.43 







(0.005E) 


(O.0517) ~ 

10 

0.33 

0.41 

1.15 

1.45 

14.7 

18.97 

20 

0.72 

0.44 

2.35 

1.47 

32.7 

21.10 

H-Sat. 30 

1 . 16 

0.47 

3.65 

1.53 

53.6 

23.67 

TiTyoming 40 

1.73 

0.53 

5.10 

1.60 

74.5 

24.00 

■bento— 50 

2.40 

0.59 

6.80 

1.71 

97.5 

25.10 

nite GO 

3.25 

. 0.67 

8.55 

1.79 

123.0 

26.43 

70 

4.21 

0.74 

10.54 

1.89 

147.0 

27.07 ■ 


*ThLe percolating solutions used were; for Na clays, UaCl ; for K clays ; 
KOI; for Ca clays, CaCl 2 ; for H clays, HCl. The normality used is 
given at the top each section. 
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TASm 7 

(After Lutz and Kemper (1959)) 

Summaiy of the Permeability averages given in table 6 and 6A 
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T A B U L K A 1 1 

Pozorovane a pocitane hodnoty rychlosti filtracniho proudeni vody v -pri 
■ h ydraullckem spadu I 

( After Kutilek 1967) 


mereno 


K kaolinit 
pocitano 

20° c 

H 

mereno 

kaolinit 30*^ 
pocitano 

c 

V 

I 

V 

I 

V 

I 

V 

I 

cm 


cm .h-1 


on .h-l 


cm. h-l , 


0 .06 

6 

0, 46 

52 

0.08 

6 

0.60 

50 

0 . 73 

26 

1,07 

101 

0. 29 

26 

1.26 

98 

0.56 

60 

1.70 

150 

0.70 

60 

2.08 

148 

0 .94 

94 

2.46 

200 

1.18 

94 

3 .90 

233 

1.41 

128 

4.28 

300 

1.83 

128 

5 . 10 

300 

2.35 

196 

8.89 

50 5 

2.90 

196 

10.27 

500 

3.37 

264 



4.35 

264 



4.85 

332 



5.79 

332 



6.33 

400 



7.52 

400 



7.86 

463 



9.31 

468 



9.35 

536 



11.38 

536 





Na kaolinit 20 

°c 


Ka kaolinit 3o°c 

0.07 

26 

0.12 

51 

0.09 

2S 

0.14 

'50 

0.14 

6o 

0.26 

99 

0.24 

60 

0.36 

100 

0.27 

94 

0.48 

151 

0.35 

194 

0.61 

!50 

0.40 

128 

0.73 

201 

0.58 

128 

0.94 

200 

0.63 

19 6 

1.39 

300 

0.77 

196 

1.76 

300 

1. 10 

264 

3.44 

502 

1.31 

264 

4.21 

500 

1.62 

33 2 



1.99 

332 



2.27 

400 



2.78 

400 



3 . 12 

468 



3.75 

468 ■ 



3 .81 

536 



4.72 

536 




Notes Hesaalte of Eutil^ (l967) as preseited in pages 75 to 76a are 
givaa in the original language of the paper( Czech.) as no 
authmtic translator was ayailable. 
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Pokracovani tabu Iky 11. 


mereno pocitano mereno pocitano 


V 

I 

V 

I 

V 

I 

V 

I 


Ca kaolinit 20°C 



Ca kaolinit lo^c 

0 .03 

9 

0.25 

50 

0.04 

9 

0 .30 

50 

0-15 

37 

0 . 58 

101 

0.22 

37 

0.68 

100 

0.46 

86 

0.96 

150 

0.59 

86 

1-16 

150 

0 .83 

13 4 

1.49 

200 

1.07 

13 4 

1.75 

201 

1.25 

133 

2.76 

307 

1.63 

183 

3.21 

301 

2.41 

280 

5.80 

495 

2.86 

280 

7.40 

502 

3 .70 

377 



4.69 

377 



5.53 

474 



6.87 

474 



7.06 

571 



9 .00 

571 




H montrnorillonit 20^ 

c 


I 

■i montrnorillonit 







3o'^c 

0 .0078 

107 

0.00286 

50 

0.0091 107 

0.003 

43 

0 .0158 

162 

0 .0071 

101 

0.017 

162 

0.0068 

86 

0.0225 

220 

0 .0125 

151 

0.0248 220 

0-0118 

130 

0.0315 

277 

0.0185 

204 

0.03 40 277 

0.0170 

170 

0.0387 

334 

0.0336 

304 

0.0438 334 

0.0311 

260 

0.0472 

390 

0.0679 

504 

0.0 550 390 

0 .0626 

430 

0.057 

447 



0 .0668 447 



0.0692 

515 



0.0793 515 




Na montrnorillonit 20^c 


Me 

a montrnorillonit 







30^0 


0.00088 

160 ■ 

0.000064 

50 

0.0023 160 

0.00045 

49 

0.0027 

245 

0.00033 

loo 

0.0052 245 

0 .00 105 

98 

0.0061 

330 

0.00083 

152 

0 .0109 330 

0.00 196 

143 

0 .0122 

415 

0.0015 

20 2 

0-0175 415 

0.0034 

195 

O .0201 

500 

0-0047 

297 

0.0 2 58 500 

0.0082 

293 

0 .0283 

58 5 

0 .0199 

499 

0.03 73 585 

0.0262 

489 

0.03 78 

670 

0.0305 

604 

0.0503 670 

0.0378 

536 
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Pokracovatii tabulky 11 (contd...) 


rnereno pocitano mereno pocitano 


V 

I 

V 

I 

V 

I 

V 

I 


Ca montmorillonit 2oS 


Ca 

montmo 

rillonit 

e 30* 

0 ,0138 

128 

0.00506 

56 

0 .0166 

128 

0.0056 

50 

0.0242 

196 

0.0117 

114 

0.0270 

196 

0.0122 

loo 

0 .03 64 

264 

0.0201 

17 4 

0.0415 

264 

0.0199 

151 

0.0488 

332 

0.0279 

222 

0.0525 

332 

0.0288 

201 

0 .0523 

400 

0,0488 

332 

0.0697 

400 

0 .0486 

301 

0.0775 

468 

0.101 

560 

0.0862 

468 

0.0963 

502 

0 .09 60 

536 



0.108 

53 6 




K illit 

20° C 


H illit 3o' 

° C 


0.0118 

107 

0.00 49 

51 

0.0146 

107 

0.0055 

51 

0 .0 20 5 

163 

0 .011 

101 

0.0228 

163 

0.0133 

10 4 

0.0274 

220 

0.0178 

151 

0.032 

220 

0 .0203 

150 

0 .0378 

227 

0 .0262 

208 

0.0433 

277 

0.0 289 

201 

0.0516 

333 

0 .043 

296 

0,0546 

333 

0.0 49 2 

300 

0.067 

390 

0.0935 

504 

0.0725 

390 

0 . 10 4 

508 

0 .080 2 

447 



0 .090 

447 




la illit 

o 

20 C 


Na : 

illit 3o°C 


0 .0056 

160 

0.0021 

72 

0.0071 

160 

0.0019 

.50 

0.0094 

245 

0.00326 

103 

0.012 

245 

0.00 408 

100 

0.0148 

330 

0.00 504 

154 

0.0173 

330 

9 .0067 

155 

0.0199 

415 

0.00745 

206 

0.0239 

415 

0.009 46 

203 

0.027 

500 

0.0131 

309 

0.0316 

500 

0.0164 

310 

0 .0336 

585 

0.028 

515 

0.0398 

585 

0.0324 

502 

0.044 

670 



0 .0518 

670 




Ca illit 

20°C 


Ca illit 30°C 


0.00605 

116 

0 .0038 

73 

0.0081 

116 

0.0031 

51 

0.0112 

178 

0.G057 

105 

0.0135 

178 

0.0075 

104 

0.0154 

240 

0.0093 

157 

0 ,0176 

240 

0.0114 

150 

0.021 

302 

0.013 5 

209 

0 .0274 

302 

0.0163 

201 

0.0 281 

3 63 

0.024 

314 

0.0338 

3 63 

0.0 277 

300 

0.036 

425 

0.0515 

323 

0 .0402 

42 5 

0.053 5 

508 

0.0438 

487 



0.056 

487 
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T A B U L K A 12 

CA LCU la T/O N of J-a 

Vyhodnoceni flltracnich krivek: Pocatf^cni bydraulicka vodivost 

hydraulicka vodivost M a Iq 


(After Kutilek 1967 ) 


Jilovy kaolinit montmorillonit illit 

miasral 


Vymenny H 

Kation 

N a 

Ca 

K 

Na 

Ca 

H 

■Na 

Ca 




1^ _ 


i_e 




-3 , -1 









M.lO cm.h 23,7 

8.1 

14.8 

1.60 

0*35 

4.2 

18.2 

0.86 

13.2 

M' . 10“-^,^ -u-l 1 c 

^ cm* n 1,3 

0.8 

1.5 

0.13 

0.07 

0.12 

1.1 

0.29 

1.5 

M/li' 15,8 

11.4 

9.9 

12.3 

5 

35 

16.5 

3 

8.8 

Iq 160 

395 

113 

171 

3 40 

29 5 

18 5 

18 5 

230 




2, 

s e r i 

e T 20‘^C 



— 2 —1 

M, 10 cm.h 27.4 

15, 1 

25,3 

0,18 

0,107 

0,26 

0*3'? 

0,10 5 

0 . 19 

M' .10~^ cm,h“lg j 

. 2,1 

4.3 

0 .05 

0.002 

0.08 

0.09 

0.025 

0.045 

M/M' 3.15 7.19 

5.88 3.6 

53.5 

3.25 

4.1 

4.2 

4.2 

Iq 200 

310 

283 

130 319 

178 

310 

29 5 : 

300 




2. 

s e r i 

e T lo'^C 



M.lO'^on.h"^ 31.5 

18.0 

29.7 

0.20 

0.130 

0.28 

0.39 

0.12 

0.22 

H'. I0“^cm.h”^ 11.0 

2.8 

5.4 

0.06 

0.006 

0.10 

0.10 5 

0.03 5 

0 .06 

M/M' 2,96 

6.43 

5,5 

3.33 

31.65 

2.8 

3.71 

3.43 

3.G7 

Iq 19 5 

300 

28 5 

120 

300 

172 

300 

29 5 

300 

^"20^^30 0.87 

0.84 

■ 0.85 

0.90 

0.82 

0.93 

0.9 5 

0.37 

0.86 

”2oA!3o 0.79 

0.75 

0.795 0.83 

'0.33 

0.80 

0.86 

0.72 

0.83 

(M/M’/M/M^O 

^q20/^o30 1,03 

1.12 

1.07 

1.08 

2.46 

1.16 

1.10 

1.22 

1,15 

1,03 

1,01 

1.08 

1.06 

1.04 

1.03 

1.00 

1.00 

^^ 2 oA ^'30 0.79 

0.58 

0.75 

0.85 

0.38 

0.83 

0.81 

0.79 

0.80 


zkmerenych 

dat 


1st Seris : I increased from O 

2nd seria : I (hydr. gradient) decreased from initial highest value 
■to zero . ' 
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state that if the entire flow cross-section is occupied by 
bonded water, seepage may not start untill an initial gradient 
is reached. The value of the initial gradient depends on the 
water content of the soil. After the initiation of movement, 
the seep^g fluid behaves like a fluid whose behaviour is between 
that of Bingham body and the generalised Newtonian fluid. 

With increasiiS^ water content, the properties resembling 
those of Bingham body disappear and the movement takes place 
in a fashion characteristic for the generalised Newtonian 
fluids. In support of this hypothesis, Karadi and Nagy (1961) 
conducted some velocity -gradient tests on clay samples of 
different water content and the results are reproduced in 
Pig. 41. Experimental results indicate that up to certain 
water content , seepage does not occur unless the hydraulic 
gradient exceeds the initial value i^ o5r-ij^_. Beyond this value, 
the relationship between seepage velocity and hydraulic 


gradient is described by a curve which joints the Bar dy^,_ line 
at a certain limit gradient (ij^). 

Initial or threshold gradient is found to increase 
appreciably with reduction in water content and the summary 
of the results are shown in Table 8. 

Swartzendruber (1962a and b) replotted some of the 
data of earlier workers e.g. Von Engelhard t and Tunn (1955), 
Hansbo (I960) and lutz and Kemper (1959) and fitted them in 
equations 14 and 15 as shovm in Pigs. 17, 18, 19 and 21. 
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Siunmary of Zaradi and Efagy’s Result (After 
Karadi and Zagy 1961) 


1 ! . t ' ! 

Exfigrlmoixt {Water {Permeability ‘Initial gradient {limit gradient 
number {content {coefficient { (i^ or i, ) { *-0 

I !(Z cm/sec.) ! ° ^ ^ 


1 

29.1 

0.72 X 

10-8 

22.8 

42 

2 

32.2 

1 .03 X 

10-8 

12.8 

31 

3 

34.5 

5.32 X 

10“® 

7.3 

22.2 

4 

39 -.3 

11.20 X 

10-Q 

3.1 

14.8 

5 

' 42 . -1 

15.20 X 

10 "^ 

0 

8.6 


The experimental results of Li (1963) for Houston clay 
(clay content 12 % and initial void ratio of 0.85) as shown 
in Pig 42 indicate a threshold gradient of around 50 and 
essentially a linesir velocity -gradient curve after that. 

Miller and Low (1963) conducted several tests on Za 
and lithium clay pastes. Their experimental graphs are 
reproduced and are shown in Pigs. 43 to 46. As seen from 
these figures, threshold gradient for Za-clay paste and 
li-clay paste is established and its magnitude being dependent 
on clay concentration in the paste. The nonlinearity of the 
flow rate-gradient relationship is not very conspicous and 
inconclusi-®B.. Miller and low state that for most materials 
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exhibiting a yield value or threshold gradient j the flow 
rate-gradient relationship is initially curvilinear (Powell 
and Eyring 1944). This curvature is attributed to any one of 
the following causes. 

(1) Reversible re-orientation of particles with the 
increase in hydraulic gradient 

(2) If there is a range of pore sizes, the yield value 
decreases in successively smaller pores as the 
hydraulic gradient increases. 

and (3) Structural breakdown of water under stress begins 
near the center of the flow channel, where the 
water is farthest from the clay surfaces and then 
progresses towards these surfaces as the gradient 
increases i.e., effective flow area gradually 
increases with gradient. 

The lack of curvature in the experimental data is explained 
by Miller and Low as due to (1) the particles being so close 
together that they could not orient along the streamlines 

or (2)' the limited interstitial volume precluded a wide range 
of yield values for the pore water or (3) there was little 
variation in the structural rigidity of the water. To clarify 
the last point, Miller and Low state that the structural 
rigidity of the water decreases with distance from the particle 
surfaces. When the water films on these surfaces are thin, 
there is less distance over which the decrease can occur than 
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when they are thick. Hence, there is a smaller total 
variation in structural rigidity of the water in the. former 
case than in the latter. The water films in the concentrated 
clays were thin relative to those in the less concentrated 
clays. 

Miller and low (1963) e.lso points out that water 
structure in clay resists disruption, but once the 'structure 
has been disrupted by flov/, it does not immediately reform 
even though the disruptive force is removed. 

To clear up the problem of filtration in the direct 
vicinity of the origin and to check the most suitable equation, 
Kutilek ( 1964 ) conducted model filtration experiments with 
minerologically homogeneous and homoionic (H, Fa and Ca 
derivatives) clay fractions and the materials used were 
Kaolinite, Montmorillonite and Illite . Kutilek chose these 
materials because of the maximum development of the surface 
phenomena in these materials, the shape of curvature in the 
velocity-gradient relationship and all the factors influencing 
filtration can be more exactly determined. Distilled water 
and the solutions of HCl, FaCl, Fa^SO^, CaCl 2 , Ca(F0^)2 in. 
different concentrations were used as the permeating fluid. 

The hydraulic gradient ranged from 0 to 500-600. The 
general nature of the experimental results as obtained by 
Kutilek ( 1964 ) is reproduced- in fig. 47. For distilled 
water, velocity-gradient relation always showed a nonlinearity 
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at low gradients whereas other solutions with high normalities 
showed essentially a linear plot. For low normality solutions 
slight curvature was noted. The experimental results were 
fitted to equation (16) 


V = 


I log (A + - I 


o 


(16) 


where 

^o 


A = (M/M' - 1) 

The parameters M, M', and are calculated from the fitted 
curves and are shown in Table 5. All the experimental curves 
were observed to pass through the origin and hence no threshold 
gradient was in sight. 

Kutilek (1966) further conducted velocity -gradient tests 
with these clay minerals which were saturated with quinolinium, 
pyridinium and calcium. The experimental results are once 
again fitted into eqn, (16). The nonline ar-ity is explained 
in terms of non Fewtonian behaviour of soil water. (Tables 
9 and 10 and Fig 48). 

Kutilek (1967) investigated the effect of temperature 
on the velocity-gradient relationship using eqn. (16). 

Complete set of this data is presented in Tables 11 and 12. 

A partial plot of these data is made and is shown in Fig 49. 
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The data shows that at all temperatures 5 the velocity -gradient 
response is nonlinear for allmost all the samples but the 
more conspicous ones are with Fa and Ca base. For all 
temperatures, the existence of threshold gradient appears to 
be certain for Montmorillonite and Illite but for kaolmnite, 
it is difficult to draw any conclusion. 

Kondon ( 1967 ) concluded that initial or threshold 
gradient increases with decreasing porosity, with increasing 
specific surface and with decreasing hydraulic radius 
(= surfa c e ^ explained these results on the basis 

of rheological properties of water in thin slits and narrow 
capillaries. 

Churayev and G-orokhov (1970) demonstrated that for unsaturated 
soils also water conductivity must be a function of not only 
moisture content but also of the pressure gradient (i,e 
velocity -gradient relationship is nonlinear). They explained 
this effect as due to visco-elastic properties of adsorbed 
water that separates the capillary wedges. 

Kutilek ( 1971 ) studied the effect of humus upon 
velocity-gradient relationship. He found he humusless 
Rendzina (loamy clay) showed nonlinearity at low gradients 
but the existence of threshold gradient was not assured. But 
Kutilek observed a linear curve for lend*ina with humus. 
Kutilek explained the results in terms of non Fewtonian water 
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JABLE 9 

Values of flow velocity V measured at hyd, grad. I . 

(After Kutilek 1966) 



Ca \ 

f 

Quinolinium | 

! 

Pyridinlum J 

■ I 

t V f 

I 

— f 
f ^ T 

f V f 

I ■ 

! T 

T V \ 

t „ t 

cm/hour 

Kaolinite 



8 

0.38 

8 

0.36 

23 

0.04 

18 

0.92 

18 

0.90 

53 

0.11 

23 

1 .49 

23 

1 .50 

131 

0.56 

28 

1 .57 

28 

■ 1.97 

210 

0.65 

56' 

3.63 . 

34 

2.33 

310 

2.72 

73 

5.02 

73 

5.55 

420 

4.55 

•95 

6.92 

129 

9.20 

495 

5.70 

1 29 

10.18 

184 

14.20 



180 

15.25 

240 

18.0 


550 26.6 





Montmorillonite 


63 

0.018 

8 

2.8 

8 

1 .5 

111 

0.032 

18 

9.2 

18 

3.6 

170. 

0.058 

28 

13.2 

23 

5.8 

506 

0.090 

40 

25.0 

28 

6.9 

257 

0.152 

56 

36.3 

33 

8.4 

302 

0.234 

73 

47.5 

56 

15.4 

400 

0.485 

129 

90 

73 

21 .6 

485 

0.84 



128 

38.8 

592 

1.28 





Illite 

60 

0.122 

8 

0.41 



136 

0.364 

18 

1 .00 



201 

0.687 

23 

1.38 



240 

0.98 

28 

1.62 



316 

1.33 

40 

2.44 



360 

1.95 

56 

3.35 





73 

4.35 





95 

6.40 





139 

9.52 





TABLE 10 


Flow parameters M (hydraulic conductance), M' (initial hydraulic 
conductance), and (intercept of M on I axis) 



(After 

f 

Kutilek 

1966) 


Slay Mineral Cation 

1 1 

t I 

* I ' 

\ 

t ! 

! t 

1 1_ 

t 

1 . 

in cm/hour 

X 10"^ 

M/M' 

1 

Kaolinite 

Ga 

113 

14.8 


1.5 

9.9 

Quinolinium 

32 

103 


50 

2.06 

Pyridinium 

4 

79 


45 

1 .76 

Montmor i lloni t e 

Ga 

295 

4.2 


0,12 

35 

Quinolinium 

14 

810 


300 

2.7 

Pyridinium 

7 

320 


160 

o 

CM 

Illite 

Ga 

230 

13.2 


1.5 

8.8 

Quinolinium 

32 

90 


50 

1.8 
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in the humusless soil and he claims that humus prevents water 
from gaining such a type of structure which makes the soil water 
to behave as non Newtonian. The experimental results are 
reproduced in Tig. 50. 

B. EXPERIMENTAL NON DARCY DATA TOR SARD AND OTHER INERT MATERIALS 

Compared to the appreciable amount of data available for non 
Darcy velocity-gradient relationships (at low gradients) for 
clays and clay minerals, there is only limited data in the 
literature for sands and other inert materials. 

One of the earliest thorough investigations with wide 
range of materials including sand is due to Slepicka (1961), He 
conducted velocity gradient tests over a large range of gradients 
and the general relationship obtained is shown in Eig. 13 and 
the fitted equation is of the form 

V = f^i^ . . . . (13a) 

Slepicka fomd the value of f as 1.3 to 2.3 in the prelinear 
regime depending on the value of K in the Darcian (i.e. linear) 
regime. The experimental variation of f in the prelinear regime 
and i- (critical gradient as explained in fig. 13) with Darcian 
iC is reproduced in Pigs. 14 and 15. 

Dudgeon (I 966 ) conducted tests on (a) giver G-ravel, (b) 
Crushed Blue Metal and (c) Glass Marbles for various grain sizes 
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and porosities. In the low gradient zone, his fitted equation 
is of the form, y = Ki^, in which n varies from 1,12 to T.27 
(Pigs. 51, 52 and 53). 

Miller, Overman and King (1970) observed non linear 
velocity -gradient relationship for fine and medium pyrex 
sintered glass frits. Phenol, added to the water destroyed 
curvilinear ity and re-established Darcy flow. 

The nonlinearity observed in sand and other inert 
material by various authors as described above is attributed 
mainly to the phenomenon of streaming potentials and electr oleine tic 
effect but they are inconclusive. 

C. M DiiRGY EESUITS ATTRIBUTED TO CAUSES OTHER THAU THE SURFACE 
PHEUOMEM 

Olsen (I960) showed the abnormality of flow behaviour in 
clay by comparing his experimental results with the Kozeny's 
equation ( 1927 ) which is based on Darcy's law. (Figs. 54 
through 57). The deviation from Darcy ' s law, i.e. the parallel 
tube model is explained by Domain Structure Hypothesis as 
explained earlier, 

Martin (I960) attributed the nonlinear velocity -gradient 
relationship, as obtained by Yon Engelhardt and Tunn (1955) to 
the clay fabric changes. Olsen (1965) holds the contamination 



dient 








95 


of capillaiy tube as responsible factor for various non Darcy 
data of different authors. 

Drom the tests on compacted silty clay, Mitchell and 
Younger (1967) observed a nonlinear velocity-gradient relationship 
upto certain gradient for both increasing and decreasing gradients 
with a hysterisis loop as shown in Dig. 58. The nonlinear 
behaviour is attributed to the particle migration rather than 
abnormal water properties due to surface phenomenon. 

Bondarenko (1968) postulated that the liquids with 
hydrogen bonding show nonlinear velocity -gradient relationship. 

He showed (Dig. 59) that liquids with hydrogen bonding e.g. H^O, 
G2H^0H, CHgOH and a mixture of CH^CO, GH^ and CHCl^ behave like 
a Bingham liquid characterised by a plastic viscosity and 
yield stress of the order of 10~ to 10“^ dynes/cm and 
hence correspondingly these liquids show a threshold gradient 
before any movement starts and once the motion is initiated, the 
velocity-gradient relationship is linear. But the liquids with 
no hydrogen bonding e.g. GCl^, OH^COCH^ and GHGl^ show a purely 
Darcian behaviour (a straight line passing through the origin), 

Childs and Tsimus (1971) contradicted Bondarenko*s 
result and produced experimental results with glass spheres 
fully obeying Dancy's, law and no threshold gradient. 

Hovak (1972) obtained non linear velocity -gradient 
relationship for Bentonite ■ and Kaolinite with a bysterisis loop 
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(Pigs, 60 and 61). FoTak attributed this to the particle 
migration and volume change during gradient application. 

D. SUMMARY OP PUBIISHEU DATA AM) ASALYSIS OP DEVIATIONS PR CM 
DARCY BEHAVIOUR 

In an excellent review Kutilek (1969) listed the possible 
errors which might have caused the reported deviations of the 
published non Darcy data. The list of errors, as suggested by 
him is shown in Table 13. 

Olsen (1965) observed that slight contamination of the 
capillary tube in the measuring device might produce considerable 
spurious data and in this light he calculated the possible 
error in some of the published data. A complete analysis of the 
reported deviations is made and is shown in Table 14. This 
table is prepared out of Kutilek’ s (1969) and Olsen's (1965) 
observations but most recent data are included by the authors. 

Swartzendiruber (1968) fitted some of the experimental 
results available into the existing two empirical equations 
proposed by Kutilek and Swart zendruber and compared the non 
Daroy effect by the ratio The slopes and 

and explained in Pig. 62. The calculated results are 
reproduced in Table 15. 
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TABLE 15 

LIST OF EXPERUffiBTAL ERRORS WHICH MIGHT CAUSE THE REPORTED 
deviations in THE PUBLISHED NON DARCY DATA. 

POSSIBLE SOURCE OF ERRORS ; 

(a) Open sample, swell with time and hence porosity changes 
with gradient application 

(h) Closed sample, swelling under constant volume causes 
gradual decrease of effective porosity with time and 
gradient application 
(o) Air trapped in the soil 

(d) Bacterial growth during prolonged test 

(e) Addition of phenol and similar solutions to remove 

bacteria might change the "quasi -crystalline" structure 
of soil water 

Cf) OantamLnation of the wall of tto oapillary tube and other 
measuring devices 

(g) Unsaturated soils-geometry of the flow path changes with 
saturation and hence gradient application. 
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TABLE l4 

NON DAUCY PWW JV PERSPECTIVE 


Investigator 


Sample 

Characteristics 


No. of Condition 

Tests of 

Per for- Plow 

mod. 


Observed Deviations 


Estimated 

by 

Olsen (1965) 


Possible 

Sourse 

of 

Error as 
Suggested 

by 

rutilelc 

(1969) 


1. 

King ( 19^8, acc. 
to Swartzendru- 

Sand and 

Sat. 

A 


. ? 


ber 1962 b) 

Sandstones 





2. 

I2bash(l931) 

Gravel 

Sat. 

A , B 


? 

3. 

Puzyrevskaya 
(1931) Acc. 
to Poliibori- 
novakochina, 

19521 

Clay 

Sat. 

A 


? 

4. 

Von Engelhardt 

Sandstbne 14 

Sat. 

B,E (?) 

0.^1^ <500 



and Tunn(l955) 

with Clay 
fraction less 
than 5 -A 





5, 

Lutz and 

Unconfined 3o 

Sat. 

B , N 

-50<lo 

a< 

b< 


Kemper (1959) 

Samples of 

Bentonite, 


E{?) 

<200 



Halloysite and 

Bladen clay 







Pastes 





6, 

Han0bo(196o) 

Natural Clays 10 

and confined 

Sat. 

B 

0<1„<2.2 

f 



sarr^le 




7. 

Slepicka(l96l) 

Sandstones 

Sat. 

B . C 


d 

8. 

Miller and 

Low (1963) 

1) Bentonite 

Pastes 

Sat. 

H , A 


f 



ii) Confined 6 

Samples of Li 

Sat . 

K M, 

N, D 

0<lt<70 




and Na Mont. 





9, 

Li <1963) 

Clay 

Sat, 

A 


£ 

10. 

Rawlins and 
Gardner( 1963) 

Silty Clay 

Loam 
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A, B 


g 

11. 
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Sif^Loam 

Un sat* 

B, E 


9 

12. 

Ha<ia8(1964) 

LoeflL Clay 

Un sat. 

B-C 



13.. 
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-and Tchuraev 
(1964) 

. Peat 
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E 



14. 

Kutilek 

Kaolinite, 
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a 
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APPMDIX B 


be read with Vol. I, Chapter 2, Sectim B) 


CAPIELAET MODEL POE SJHPACE ACTIVE SOILS 
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Yelocity distribution, : 

Por the itee mter zone, Integration of equation (2.47) 
yields 


~ 4 °2 


(2.67) 


where 


Cg = constant of integration 

For the bound water zone, integration of equiion (2.59) 


yields 


* " 


e [-1 + I, (A ) T (X) 

*■2 1 1 o *^ 4 


% 


j- I^(^)l+ Cj 

(2.68) 


vdaere 


Cj = constant of integration 


12^^) > / X e‘^^dX = - 

( A ) 

->-rr- i>. -V [(GX-1) logx 


1 + AX -AX 


A 


(2.69) 


GA-| 1 GX 


X e 


AX 


BG 


G 


G = B - A 


2 o 

V?- r 

V 3 "f 

B 


2 '• i . n!n •* 
n=1 


( 2 . 70 ) 

(2.71) 

(2.72) 

(2.73) 
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I^(x) and are defined earlier In. equatims (2,52) and 

(2.54) 




= logX+ I 
n=1 


^ equations (2.39) and (2.68), one gets 


Cj = e 


AX I ( X ) i 

o- r 2'' s"^ w 






where 


(2.74) 


(2.75) 


^ 2 (^ 3 ) * value of I (X) at X = X 

s 

^(Xg), » Talue of I^( X) at X=* X^ 
^^(X ) = Talue of I (X) at Xb X 

*r « 4 g 



(2.76) 

(2.77) 

( 2 . 78 ) 

(2.79) 


l 2 ( Xg ) » Talue of I^CX ) at X = X^ 
^ = Talue of I 2 (X ) at X a X 

•mm ^ O 

t . 

^4 (Xo) - value of (X ) at X = X^ 


(2.80) 

(2. 81) 
(2^82) 
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Equations (2.6?) and (2.68) along with equations (2.75) and 
(2.79) complete the solution for velocity distribution. 

Average velocities 

Average velocity over the entire cross section of the 
capillaiy is obtained from eqn. (2.83) 


1 r f 
TTR O 


u 


E 


2'trr dr + / u 2irr dr ] 

r ^ 

0 


Equation (2.83) can also be written as 

♦ = 2 [ / «_XclX+ / 4) XdX ] 

0 X 

o 


where 


(2.83) 


(2.84) 


u 


^av 


av 

Th^ 

0 

Vo 


(2.85) 


aibstitutlng expressions for 4^ and 0^ from equations (2.67) 
and (2.68) in eqn. (2,83) and integrating one gets, 


d) = 2 [ 


^2 ^ i _ ^2 ^ 

16 2 o 2 s 0 ^ 


O 


(2.86 ) 
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Equation ( 2 . 86 ) can be 


u = i ( J. 0 

av '■ Q 


rearranged as 

2 ,4 

c 

AX 


YqE 

8 ~ ) [ (BOjX^ - ) + 8 { o,(xJ - x2 ) 


3 s o 

“(l 5 (X 3 ) -L(XJ )>1 


f Y H A 1 

2 J ( l£_ ) 


Y 


o 1^0 


-<v\)- V^o)” 


( 2 . 87 ) 


where 


defined by equations (2.75) and (2.79) 
Ag is related wilh i by equation (2.62) 

l 5 (x) = / X 12(A) dX 


^ 2 2 

~ 4~ [ A X + 3AX + 3 ] 

A 

33 (Xg)= value Of I_(X) at x= X 

^ s 


(2.88) 

( 2 . 89 ) 


^5(Xq)= value of iJX) at X=X 

^ o 

= / xi^(x) ax 

= f- [logX-i]+ x^ f (_1)“ - (Axf 

n=1 n(n+2)|n 


(2.90) 


(2.91) 


Ig(Xg) = value Of I (^) at x= X 

^ a 

at X= X 

O o 


(2.92) 

(2.93) 
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(x) = / X I^(x) dX 


a X 

(ff - 20). + 2) 


- ( ) ( -^ e®’' ) 

• nrx ^ VJ 


+ -— [ log X - -1 ] 


+ p x^ [ y -v-- ^^/., - ] 

^2 n(n+2)ln ^ 


(2.94) 


I^(Xg) = viue of I^(x) at X =5 X^ 


(2.95) 


I^( X^) = value of I^( X) at X = X^ 


(2.96) 


and Pg are given by equations (2.72) and (2.73). 

Equaticaa (2.87) completely defines the average velocily 
gradiait respmse for the system. It is evident from ecji. (2.87) 
that velocity gradient response remains non linear as Icaig as 
Xg< 1.0 or i < i^. Once i exceeds i^^^ the flow remains 
fully developed along the entire cros section of the capillaiy 
and the V elociiy gradient response takes a linear shape of the 
form 


1 '^o® 

^av 8 ^0 ^ ^ ^ ^ ^0 ^ ^ 


(2.87a) 


where a and b a3?e functicais of m md X , 

0 
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In the absence of the bnund water zone 


X 

0 

= X = 1 
s 

(2.97) 

f 

w 

= 0 

( 2 . 98 ) 

and 

A = CC 

(2.99) 

Hence I (x ) = 0 

ii S 

(2.100) 

03 = 

0 

(2.101) 

°2 = 

1 

T 

(2.102) 



(2.105) 



Ill 


AA 


P. = 


^2 


"3 = 


2e 


AA 


[ (t+A) 


8 


•“AX 

(l+AA ) e °] 


2 £ 
A^ 


AA 


2e 


~A, 


A 


2 ^ [ (1+A) e "j 


(2.106) 


(2.107) 




and as foiind earlier 


(2. 108) 


A 


log m 
1-A 


(2.31) 


Equation (2. 104 ) completely defines the velocily gradient 
response in the oapiliaiy element with bound and free water zone 
having an exponmtial variation of viscosity within the bou:.d 
water zone but without the presence of any interacting force. 

If the soil IS assumed to be composed of bundle of strai^t 
parallel capillaries of untfoim radius R and using the concept 
of iydraulic radius and equivalent grain diameter 
Dg(!feylor 1948), equation (2.IO4) can be written in the foim 


V ' t 72 “s if? V 

o ^ 


(2,109) 


where e = void ratio 

!Ehe factor sqp.. (2.109) has come because of circular 

capillaries and hsice using a general shape factor 0 for an 
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irregular capillaiy, eqn. ( 2 . 109 ) becomes 




( 2 . 110 ) 


Oomparing this with Darcy's law, one gets, the expression for 
coefficient of permeability K as 


K 


= [ G D 
'■ s s 


1 +e 


] { 4 P(m, X^ ) } 


(2.111) 


Equation (2. Ill) gives the permeability or the hyirailic 
conductivity coefficient taking into ccaasideration caastant 
viscosity in the free water zcaae and exponential increase of 


viscosity in the bcund water zone of a soil water system. But 


the terms in square bracket in equatioi (2. Ill) give Kozeny's 

hydraulic conductivity coefficient (K ) for a soil water 

kozeny’ 

system having ccnstant viscosity throu^out, theorefore, equation 
(2.111) can be rewrittm as 


K 




(2.112) 


Lozoay 


The ratio g- may be looked on as a discrepancy ratio 

nozeny 

arising out of the existaice of two distinct zcnes of water 
(free water zone and bound water zcne) and the exponaatial 
increase of viscosity in the bound water zone. 

As a special case, if free water zone is absent (i.e, Xq= O) 
in a clay water system, equation (2. 112 ) reduces to 
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ozeny 


= 44'(ia) 

no fjpee water 


where 


(2.115) 


= P(m,A^) (2.114) 

Xo = 0 

The function ^(m) is obtained by putting X^ = 0 in the 
expressicaa for ^(m, X^) (equation (2.105 )). 

Equations (2. 112) and (2 . Ill) clearly bring out the effect 

of chan ^d visoosily of the bound water on the permeabilily 

coefficient K and they are plotted in figs. 2.7 and 2.8. 

Pig. 2.7 shows the variation of the discrepancy ratio (ratio of 

actual pemeability K and the permeability given by Kozeny ’s 

eqn. ) with the change in the values of m and X^ and Pig. 2.8 

shows the variation of discripancy ratio with the change in m 

values for a system where all the watei px'esent is only bound 

water. is seoa from the derived equations end their plots 

that even within the Darcian linear zone of velocity-gradient 

response, the hydraulic conductivity characteristics for surface 

active soils, among other things are also a functions of m . aid 

X^, The actual permeability is found to be veay much sensitive 

to X but some vdiat less sensitive to the m values. In the 
o 

absence of bound water (Xq=’ 1) (as may be expected for sand aid 
other inert material) the discrepancy ratio becomes 1.0 as 
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expected. The practical ran^ of ’m’ and may vary from 

1 to 500 and 1 to 0 respectively depending on the clay content 
and porosily of the soil. 

In the absence of any free water zone (as may be expected 
in case of comparitively dense clays of high surface activily), 
the discrepancy ratio is observed to decrease very rapidly upto 
’m’ values of 50 and the rate of decrease becomes smaller 
there after. 
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APPMDIX - G 

WATER MOVMEKT CAPILLARY ZdTE 

Movement of saturatim line in the capillaiy zone is of 
interest in earth dams and other similar structures* Assuming the 
flow to be one dimensional (either hoiizoital or vertical), a 
theoretical expression for the movement of the saturation line as a 
function of time would be of great help to the field and desi^ 
engineers for estimating the time dependent strength and stabilily 
of the capillary zone. When the flow is assumed to be I^rcian, 
expressions for both horizontal and vertical capillary rises are 
available (Taylor, 1948). Here, an extension of the same for ttie 
non Darcy flow defined by equation (C.l) is made and presented. 
Bxpressiaas for both the rate of horizontal and vertical 

V = Mi*^ . (c.l) 

Capillary saturatiaa are forwarded. 

RATE OE HORIZCITTAL CAPILLAEZ SATCRATICXf 

To assure a truly one dimensional flow, a system as shovm 
in fig. 0.1 is considered. The system, as shorn, is used many times 
for the capillarity - permeability test (Taylor 1948). Here, the 
soil sample has been brou^t to state of dry powder and thoa packed 
into a glass tube having a screen over one end and, a voated stopper 
at the other end. The tube is then immersed in shallow depth of 
water in a horizontal position and the water proceeds into the 
interior of the soil by capillary actim. The distance of saturation 



118 



Pig. Cl; SET UP POR HORIZONTAL CAPILLiJlY SATURATION 
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X as a functicp of ticie is to be found out for a non Darcy flow law 
of the type as given in equation (C.l). 


Dor the systeia as shown in fig. C.l and considering the flow 
as given in equation C.1, equation 0.2 can be written as 


S EL V 
e B 

where 


h + h ^ 
^0 0 


( 0 . 2 ) 


S = degree of saturation 
n = porosily 

G 

V = Seq)age velocily 
s 

h = Capillaiy head vihlch is constant for a given soil at a 
given void ratio 

h = As shovsn in fig. C.l 
0 


Eqn. (c.2) can further be written as 




(0.5) 


Separating the variables and integrating, one gets 


Sn 1 

e 


n+1 


(0.4) 


where C.j is the constant of integration. 

The initial condition is given by equation (C.5) 


t = 0, X = 0 

By equation (c.4) and (0.5) 


(0.5) 


c^ = 0 
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Hence, equaticaa (C.4) can be written as 

n+1 


(C.6) 


where 


Tg = Tine factor = 


Mt 
H Sq 

e 


and X = “ 
H 

H = h + h 
o c 


(C.7) 

(G.8) 

(0.9) 


Equation (g. 6) gives the desired relatiaa between time and dist^ce 
of capillary saturation in a non dimoasional form. Eor Darcy flow 
(n = l), equation (c.?) reduces to the known parabolic 3 ?elaticai 
between X and "Og as given by Taylor (l948). For various values of 
'n‘ equation (C.6 ) is plotted for X ranging from 0 to 1 and is shown 
in fig. 0.2. For a laboratory set up and for a normal size of the 
sample, the value of X will hardly exceed uniiy. It is concluded 
from equation 0.6 its plot that compared to Darcy flow, nen- 

Darcy flow leads to slower rates of capillary saturations if n < 1 , 
and vice versa. 


RATE OF YERTICAI CAFILIARX SAIUEATICW 


For this case, the ^stem along with heads acting cn it is 
schematically presented in fig. 0.3. 

At aty Instant of time, if z is the hei^t of saturation 
line, thoa by the flow law assumed (eqn. C. l) 
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. h -2 

a, ^ = 

e dt z ^ 


The abOTe equation can be non dimensicnalised as 


dZ ^ (1“Z\n 

dT ^ Z ^ 
r 


where 



_ Mt 

'^y h Sa 
c e 


(G. 10) 


(c.ll) 


(g.12) 

(C.13) 


As seen from fig. G.3 the value of Z will not exceed unity in 
this case. Separating the variables from eqn.(G.1l)and integrating 
and using the initial conditim = 0, Z = 0, cne gets 



^ n(nH-1 )(n+2)(n-f3) A 
* (n+5)Ii; ^ ^ 


(G.14) 


Bor non Darcy flow equation (0.14) is the required relation 
between time and capillaiy saturation hei^t in non dim^sional 


form. Equation (c.14) is plotted for various values of non Darcy 
exponent ’n' and shown in fig. C.4« When n = 1, equation (0.14) can 
be writtai as 



or 


» V 



(G.15) 


= -Z - log (1~Z) 
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FIG. 03: DEFIITIIIOI'T SKETCH FOR VERTICilL CAPILLARY SATURATION 
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which is the well knowQ solutioa for Darcy flow (lay lor 1948 )• 
Equations (0t14) and (c. 15) and their plot (fig* 0*4) suggest that 
for n <1, compared to Darcy flow, non Darcy flow predicts slower 
rate of capillaiy saturation at early times and faster rates of 
capillaiy saturation at later times* ffiie trend for n > 1 is just 
the opposite* 



